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ABSTRACT 

This  thesis  is  an  extension  of  the  more  comprehensive 
study  made  by  Kenneth  V/.  Goff  (5)   of  the  applicability  of 
correlation  techniques  to  the  field  of  acoustical  measure- 
ments. 

The  ability  of  correlation  techniques  to  separate  an 
acoustical  signal  at  a  point  into  compotnents  according  to 
(1)   the  source    (2)   the  transit  time  frcm  the  source  to  the 
point  In  question  ajid   (3)   "tiie  frequency;  vas  used  to 
measure  the  transient  response  of  a  reverberant  room.     The 
results  of  this  correlation  analysis  were  compared  v.-lth 
pulse  measurement  data  talaen  under  similar  experimental 
conditions.     The   results  of  the  two  methods   v/ere  found  to 
agree. 

A  techniqiie  for  measuring  the  directive  properties  of 
the  V7ave  fronts  associated  v/lth  soiond  fields  v;as  devised 
and  tested.     It  i^ras   shown  that  the  property  measured  cor- 
responded ^.^th  the  definition  of  "diffusion,*'     Results  of 
diffusion  measurements  talaen  in  a  reverberant  room  vrltb. 
various  portions  of  the  v;all  covered  with  a  highly 
absorbent  Fiberglas  curtain  were  presented  together  with 
corresponding  transient  responae  data. 

The  above  mentioned  experimental  measurements  v/ere 
parallGlod  by  an  explanation  of  ihe  theory  and  a  sample 
of  the  mathematics  involved  in  the  correlation  analysis 
of  bounded  sound  fields. 
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1   INTnorJUCTIOK 
1,1  HISTORY 

Tho  word  "oorrftlatlon"  Innlloo  tlio  act  of  ooeklns 
rslatlonnhlpo  betvroon  two  dlccroet  quantltleo  or  functloiui 
of  time  or  moro  properly,  two  oncembloo  of  ovento.  Th© 
concept  Id  oloaontary.  For  that  roaoon.  It  lo  la::>ooDlbl« 
to  pln-polnt  ooia©  exaraplA  and  elto  It  ac  "belnc  tho  flrot 
Instanoo  of  tho  applloatlon  of  correlation.  It  Id  rruf- 
floiont  to  oay  that  the  proooss  la  one  vrtaich  owec  ItB 
dovQlopMont  to  nian*0  innato  curloclty  coupled  vrlth  his 
ability  to  reaiaon, 

IHctorically,  correlation  io  a  very  important  pro- 
oeno.  ./hen  man  beoanie  dissatisfied  vrlth  the  more 
rocordlnc  of  ovonto  and  first  attempted  to  explain  then 
in  to  ma  of  a  preceding  series  of  incidents,  he  \mB   re- 
eortlns  to  a  form  of  corrolation, 

Corrolatlon,  lllw  many  other  proceooeG  involving 
basic  concopto,  lends  iteolf  to  coneice  and  explicit 
expronolon  by  mathematical  symbols.  "The  tranBlatin^  of 
the  prooeoG  from  v;ordo  to  mathematical  symbols  moro  or 
laea  parallels  the  evolution  of  the  prooeos  from  on  art 
into  a  science . 

tSie  nathenatlcal  exprensicn  of  tho  corrolatlon  pro- 
Mae  xma  a  natural  outgro^fth  of  the  dovolopmont  of  the 
field  of  statistical  methods,  Before  160C,  no 
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aathma'tloal  oonoeptiasio  of  probability  yoro  rooosnized. 
CKillford,  In  hln  book  on  poyclioootrlc  aethods*  Bays, 
•daablero  had  omculp.ted  much  conflomlng  c^anee  of  chanoo 
iftion  It  cama  tlsia  to  oonsider  thoir  Ioooog  and  Gains « 
partlcoilarly  tholr  loooos.     They  ov©n  attoaptod  to  Intor- 
est  nathooatlclons  In  thoir  problems,  thotach  vdth  rjaall 
oucc©88«     ...     Ttxe  seventh  century  saw  the  beclnnlnc  of 
eerloufl  intorost  In  the  aatheimtics  of  dhance. 

•Bomoulll  (l63^f-1705)  publlcihed  tho  firct  book  to 
\m  entirely  deyotod  to  tho  otibject*     DeMoivre  (1667-1754) 
my  b©  crodltod  ;/ith  the  discovery  of  tho  noroal  dlntrlbu- 
tlcn  curve  at  about  1733.     ^an  that  tlino  on,  Intoreot  was 
arouodd  amonc  antronooers  as  well  as  mathesiatlolans  •     By 
1812  I^plaee   (1749-1027)  had  vnriLtten  v2iat  la  conrldored 
the  rroatect  oincile  ^?orl:  on  probability.     In  It  he  cave 
proof  of  tho  Tijcthod  of  lORot  cq-uaros. 

'•It  was  GauDs  (1777-1854)  who  denonntratod  tho  croat 
practical  valts©  of  the  normal  curve,  dhowinc  how  it  ap- 
plied to  the  dlotributlon  of  lasastirsBnnts  azid  to  orrors 
ttade  in  aelentlfic  obc^srvatl one .     it  was  h©  ^o  devioed 
the  f«*3aroontal  xaathode  of  ocnputatlceti  of  aeans,  probable 
errors,   and  the  111^.     ... 

•B»  applicr>.tlfm  of  the  noroal  eiirVD  and  elementary 
statistical  Ejethodc  to  bioloGlcal  and  social  data  nust  be 

*!teforence  ntnbors  rofor  to  tlio  Dlbliocraphy. 


attributed  flrct  to  Quetelot  (1796-1079),  royBl  aotronomer 
to  tii0  JcluG  of  Bolglum.      . . .     ^Ir  PTEunois  Oaltcai  (1822- 
1911)  In  uorldJiG  on  tho  probloao  of  human  liarodlty  founfl 
that  the  nomal  cuxire  euid  Ita  siatpler  appllc-^.tlono  wer» 
inad^quatd.     He  invonted  a  nindser  of  additiooal  statietieal 
toola,  among  th&m  tha  nethod  of  carrolcition.** 

In  IBV.B9  olovon  years  aftor  Galton'e  flrot  publica- 
tion of  the  ecnoapt  d^wlopad  In  tamo  of  tho  Idaae  of 
resronclcai  linos,  tho  term  **corrol?.tlon*  flrnt  appeared  In 
print.     Pearcson,  Bdaewortli,  and   ieldon  furtlior  dovelopod 
the  aathomatieo  of  correlation  to  a  point  v;hcro  it  ifaa 
©anei^lly  adopted  by  iho  otatistiolGiis  in  tho  fioldo  of 
••OttoialcB,   nociolOGy*  biology,  aatranony,  meteorology  and 
a  few  othor  phyoioal  oolenoea*     field,  in  hie  bool:  on 
otatiotlca     rovi©i/s  tho  variottB  appllcatlona  of  the  theory 
of  oorroliition.     As  an  addendum  to  his  treatment  of  tho 
a«i^Joot«  Held  compilod  a  bibliography  \^ileh  io  ©xoollent 
in  its  Goope  of  tho  oubjoct  of  eorj^latlon. 

BoTtia'DO  the  most  Ifflsortnnt  admmfie  in  the  flold  of 
ataticticc  haa  boon  tho  natiiciia.tlcal  expositicn  of  the 
3?elatloairJiipo  betv^^on  tho  flold  of  oocimunl cation  onslnoor- 
3U3i6  and  otatlr.tic?.!  treatment  of  tii!»  ocrlon,     Tho  rigid 
mathematical  proofs  of  liorbcrt  I'^ioner^  not  only  put  the 
prooess  of  correlation  on  a  fixm  math<^atioal  foun3cition, 
but  alco  tended  to  accolorato  further  advancomontn  in  the 
field  of  oosmuni cation  engineeriiiG . 
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1.2     TIME  SERIES  AND  CORRECTION 

MtonsT  vioflnoc  tlsM  eerios  as  •?3eqtwnoes,  dlrcrot© 
or  cc3(ntlnuouo,   of  quantitative  data  aasiGned  to  speoifie 
mmmitn  in  tin©  and  otudiod  vdth  reopoct  to  tho  ntatiotloa 
of  tlielr  distribution  in  tiiao.     They  nay  bo  olmplo,  in 
i^ioh  caao  t^ioy  canaiot  of  a  single  nunerie&Uy  given 
observation  at  ea^  mccMnt  of  tho  diocrote  or  continuoue 
baet  oenuonco;  or  nriltiplo,  in  whleli  ease  they  ctmsist  of 
a  nisnbor  of  aoparat©  quantitloo  tabulated  accordinc  to  a 
tiae  eoamon  to  all«* 

The  apnliof^.tion  of  rtatictlcal  motliodo  to  cooatunica- 
tlon    ensincorlnG  roquiros  a  model  of  an  explicit  functicai 
of  time  tihooe  properties  can  b©  deocribod  in  torroo  of  itai 

1.  mean  valiie 

2.  aean-oquape  valt© 

3,  poiwer  spectrum 

4,  correlation  fimctic»i 

The  function  of  tlno,  bo  sides  being  capable  of  deccription 
in  tenas  of  a  suf fieiently  complicated  set  of  probability 
dKisities,  mtBt  also  be  a  etatlcnary  proeeosi  that  is,  one 
in  uhlch  none  of  tlie  probability  dictributlcfin  viiich 
deacribe  the  prococo  drm^  vdth  tliae. 

It  is  neeeseary  at  tliie  point  to  define  tho  ntatie- 
tleal  quantities  aecociatod  vdth  the  aeamxreiaont  of  randcn 
functions  €£  tluis  cucfe  as  the  thermal  voltage  of  vacuum 
tubes. 


■II 


The  nean  value:     If     f  (V)  dV      Ic  the  probability 
thit       V(t)       lios  boti/eon       V      and       V  +  dV      then  the 
aean  valte  of  a  noloe  Is 


<v(t)>      =         jv  fW  dV 


vnder  the  restriction  that 


•o 


f=»(\r)   dV  a  1 


■'-oft 


The  aean-sqtiar©  value  Is 


^o6 


<(7»    (t)>  =  ^'^  f  (V)   dV 


•_  oO 


^ere  the  notation    <      ]>  means  tlio  otatlotloal  averaco, 
as  comiiuted  from  the  probability  denoltloo  vfixXch  deccrlb© 
the  onoenble. 

The  time -averaged  mean  value  of  a  no!  bo  voltaige 
y(t)       Is 


vH:)  *  11m      2w  r  V{t)  dt 


The  iQoan^square  time -a  vq  raged  value  Is 


Wit)  -  llm     1,  V«   (t)  dt 


T-^*    2T 


Cr 


A  rnjadom  function  may  be  assnaaed  to  vanish  outolde 
of  the  tlae  Interval       -T/2  ^  t  ^  T/2      In  v*ilch  c^,Ge 


V(t)  =  Tuf)     o     ^^^"*^    df 

-/-Oft 


I 


^ihere   A(f )   Is  a  oonplex  function  and  •qualo  the 
volta^  spectrum  «f  V{t)         .  >  Itli  the  old  of 
Paroeval'o  theorem: 


-        '  X 


t)dt  = 


7-{t)dt  = 


/    (f )/  •  df 


''-cO  -^-yz.  ''"•* 

we  may  obtain  an  expreaclon  for  tho  total  (finito)  power 
of  a  zTindom  function: 


■7/z  r^ 


lim   JL       (    V'»(t)clt  =  "FTET  =     I '(f)  df 
Whore      '■(f)       ,  til©  poR/er  Gpootna,  is  defined  as 


11m    2  /A{f)/^ 

*Rie  auto- correlation  fwictlon  (P,,i'^)  of  a  random 
function  io 

cP„W=;  v(t)  v(t-r)  =  <^v(t)  V  (t-r)> 

similarly,  the  cronc-corrGlation  function  of  tv;o 
reaadom  fimotlono  is 

^z'^'^^    =:     Vi    (t)   Vo    (t-r)      =    <V3,    (t)   Va    (t-7')> 

Tbe  link  bettfoen  the  correlaticm  function  ar^i  tho 
other  statistical  quantities  lined  to  descrlbo  random  fitnc- 
ticjn  has  Ita  basis  in  ..'iener*G  tliooi^a:     tlie  aiito  correla- 

The  Gonoral  rolatlcaiships  between  the  function  of 


tine,  the  opeotral  denalty,  the  pamv  speetruQ,  and  th« 
oorrolatlon  function  ap«  rtiown  In  a  ctdicnatlc  devinod  by 
atnneth  i^?.  ooff^,  Fis-  i.l. 

Tho  rodulta  obtained  by  tlxo  cori'^XatlCRi  analyola  of 
acouctloal  problorao  ax^   similar  in  corao  ronipeota  to  rocrultn 
lllitalnod  by  both  ntoady  ntato  nothods  and  by  puloo  nothoda, 
Frcn  a  rudlmontary  imowlodGe  of  the  ziattiro  of  sijieedi  oaaA 
muolcal  GoundSp  on*  seta  tlio  Intultlvo  foollnc  thr.t,   olnee 
the  coundo  noi:^3ally  oncountored  in  aoountlcal  problona 
roGorablo  pulned  \ia.ye  trains,   otreosins  Uio  pulno-lilse 
aepeets  of  correlation  raothoda  woiad  ylold  aisnlflcant 
results,     IfO'w'ovor,   ouch  a  oupponitlon  lo  not  sufficient 
to  Jiictify  the  noGloct  of  the  steady  otnto  .anpoctn  of  the 
corrolr.tion  proooos.     Qie  ability  to  apply  a  docroo  of 
fr^uoncy  flXterinc  in  t2io  couroo  of  coQinitlnc  tlio  cor- 
Vilatlc^  fun0tl<m  has  aevoral  8u!vBnta@e6*     Paraaouat  opting 
theoe  covoral  advrsntageo  la  tlw  fact  that  nuch  of  tho 
jprescnt  body  of  laaoi^ds®  cai  ftttoitttieal  probloniG  Ig  based 
Oft  toots  aade  with  steady  etato,  frequonoy  filtered 
noloeo,     !iMs  fact  smd  tho  connoquent  advantaeeo  Gained 
by  the  ability  to  eenoparo  rooultn  \rii^  thoce  obtalnod  by 
more  clacclcal  oothede  maJm  it  deairablc  to  oomprcoico 
the  tlsie  flltorins  properties  of  eorrelatlon  methodo. 
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1.3    TIKE  FiL?ERina  73  PF®'        :  mv     . 

Many  of  tho  applloatloain  of  oorrolatlon  toohnlqtxoa 
In  tho  field  of  acountloa  depend  upon  tha  ability  of  th© 
method  to  ceparate  Gi|;nalo  Into  ooaponsnts  aocordlnc  to 
thoir  trrjinlt  tliao  fro©  a  so\ar<^. 

For  Inctojico,   rmpposo  tt^o  nlorophonoo  v;©ro  placed  In 
tbo  oouod  field  of  a  point  source  radlatlnc  randoQ  nolo© 
Into  fr©0  omo©.     If  tho  nlcroi^hotno  outimtc  vToro  nil::©d, 
wad  this  c<»ablnod  voltaeid  mjgvo  crooc-corr*Glatod  with  th0 
loudcpealsor  input  voltogG,  the  follotNfinG  ©ffoate  would  "be 
noticed!     At  A  dtelcty  tlnw  oorreapondlns  to  tho  trnnolt 
iteft  frcn  the  loudspeaJcer  to  the  microphcan©  cl orient  to 
the  r>o\iro©,  tha  earrelatlosi  function  wo^d  oxMblt  a  ix>ak 
value.     At  a  do  lay  tlao  oorreflpcndinG  to  the  trannlt  tl!!K> 
from  tho  lotdcpeal^r  to  tho  cocond  tnlcrophono,  another 
pMk  wcKUld  appear  In  liio  plot  of  oorrelrtl on  function  wb 
tisto  delay. 

If  the  micror^honoa  vsero  placed  too  does©  tocetlxcr, 
the  rooultant  plot  of  eorx^laticn  function  vn  tinic  ("'Gl^.y 
would  be  dlotortod.     Tho  aiaallest  difforenoo  in  tmnnlt 
tlED  -./liich  produces  an  tmdlctorted  nain  p©ak  in  the  cor- 
relation  fuacticai  correo^xxuls  to  the  reDolvino  xxxjor  of 
the  cyotcn,     QSiia  ainimsa  transit  time  siay  bo  of  tli©  order 
of  on©  period  of  th©  avomotrlc  aosoi  frequency  of  the  band« 
Tddtli  of  the  random  noln«  uc;©d|  snd  the  oorreoponding 
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reoolvlnc  pow»r  of  the   syotom  lo  equal  to  the  product  of 
thlo  nlnlmum  trannlt  tim©  and  the  cpeod  of  sound  In  tiao 
air. 

Since  the  auto-correlatlan  ftaictlcn  and  the  poiier 
STXjctrum  of  a  stationary  random  tinio  function  aro  Fourlor 
trancfom  palro,   tlioy  exhibit  certain  charactorlntic  in- 
verse nproadtnc  proiiortlea.     For  exanplo,   a  froquonoy 
spDctrum,       F(g>)       ,   flat  from       0  to  i  rx/r    has  as  ito 
trancform,   a  function  of  time,     f(t)      ,  havlns  tho  genoral 

shax^e  of    1  nin  $- 1         v;hlch  has  on©  pi-^doolnant  peak  and 

r 
axia  crossings  at'-'i^,-2T^,i4T,    ...  ,   Plg.  1,2 

(a).     A3  tho  bandi/idth  of    F(cd)     Incroasoa,  tho  axiG 

crossings  of     f  (t)     occur  at  progreoolvoly  smaller  valijos 

of  tlEi©,     Tho  limiting  cao©  lo  an    P(m)     equal  to  1.     The 

Fourlor  tranoform  of  such  a  spGCtrun  is  the  co  called 

•Iraptilne**  or  "de//a  function"  which  1g  a  function  of  tiae 

havinc  an  Infinitely  craall  duration,  an  Infinit©  ampli- 

tud©  and  an  area  of  1  under  tho  curve;  Fig,  l,2(b).     Ccn- 

voreely,  as  fhe  bandwidth  decreaoes,   the  poalis  in  the 

fimctl  on       sin  ^        tend  to  equs.lise.     In  the  llmitins 

X 

cane,  F(c})  is  coQpoaed  of  on©  frequency   (-  qo)   , 
and  the  invero©  Fourier  trraiaform  f  (t)  ie 

coDlne    (cDo)  t    J  Fls»  l#2{e). 


ill 
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Figure  1.2(a)   Fourier  transform  pair  for  finite  band- 

width  noise. 
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Figure  1.2(b)   Fourier  transform  pair  for  infinitely  4 

wide  bandwidth  noise. 
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J'igure   1.2(c)     Fourier  transform  pair  for  a   pure  tone. 
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TpmL  Hi*  aboT'e  cUseuosidi  tnd  lllUGtratlan  it  Is 
obvious  tliat  tbB  pre roqul sites  for  lii^ih  resolutlcn  in 
tin©  are  diametrically  opponed  to  the  reqiiirtjments  for 
hiGh  frequency  reoolution.     It  is  eqiJially  clear  that  most 
acountical  meacuromonts  made  by  correlation  notiiods  ro- 
qulro  a  comprcjcilao  between  th©  frequency  and  time 
filtorlng  properties  of  the  process. 

For  Inotanco,  in  measurins  transmiGoian  loos  of 
panels,   one  generally  decirea  to  Imow  ho^:  tliic  property 
varloa  i-dth  froqi©ncy.     3iufl  the  bandvddtii  of  the  tcot 
signal  auGt  be  cut  to  a  degree  cocLmennurate  vdtli  tlio  fre- 
quency roGolution  desired,     This  increao©  in  froqiioncy 
resolution  munt  not  be  so  great  as  to  create  confusion 
betv/oon  the  ouccensiv©  peaks  caused  by  the  signal  trans- 
mitted directly  throuj^  tli©  panel  and  -Uiooe  peaks  caused 
by  sicyialc  arrivlns  at  the  test  point  by  different  paths. 


i 


i 


13 


II        INnTROMffifTATICil 
2,1     INTRODUCTICH 

Baa  BUbjoct  of  InntrucKsntation  mxict  Incluao  can- 
oidsratian  of  tho  methodo  of   (a)   eathorinc  data,    (b) 
roduolng  data,    (o)   preoentlnt:;  and  rooording  resulto. 

Squimeni  for  G^tJiorijaG  data  may  Include  Bubh  itemo 
gys  electronic  nolGO  voltage  g^s^nerators,   1  oud  ope  alia  zro, 
mlcroimonea,  and  in  a  case  v/here  the  reduoins  eqixipoent 
la  not  portable,   rocorci.inG  instnsaents,     A  reducinc  systom 
is  generally  conpiiaed  of  a  nmber  of  cepiarate  electronic 
units  eada  of  t^iidi  performs  a  apocific  process.      :5ivialon 
of  the  reducing  proceGs  into  a  set  of  sub-pi?oooco03  lends 
versatility  to  the  oystea  and  also  facilitate o  testinc  of 
th©  ooaponentG.     The  end  product  of  tho  roducticn  procens 
aay  bo  proconted  in  varlouo  vrayo*     If  the   roductlcttx  nyatem 
is  electronic,   tlio  resulto  laay  be  preaented  by  meano  of 
Rjetor  road3.n£r.  or  OGClllOBCope  tracec. 

It  a  reducing  eyetem  is  alreadj^  in  oxlctenco,  tlicn  the 
equipBaent  for  Catherine  data  and  presentlns  reeultn  must 
be  selected  in  crach  a  way  ^^at  the  beot  use  will  b©  made 
of  the  features  built  into  tho  roducins  systen. 

2*2     Tim  AKAI.OG  COHEELATOR 

As  implied  by  th.e  er'U^?,tion  defining  tlio  croso- 
oorrelation  function 

<^  21       (r)  =  lim  1        \  fx  (t)  fs  (t-7^)  dt| 
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*  device  cr,  )ablo  of  connitlng  l^ils  function  rauct  porf oirm 
tlae  oporatlans  of   (a)  do  lay,    (b)  multiplication,   and   (c) 
IntGcratlon.     If  the  fxmctlons  of  time  w»r*  oloctrlcal 
voltaGos  ouch  as  thoce  produced  by  &  olcrophono  In  a 
oaund  flold,  the  application  of  thece  voltageo  dlrootly 
to  the  input  terminals  of  an  olectrlcal  analog  computer 
AooiC^nGd  to  porform  th©  roqulrod  operation  would  ylold 
tiio  corrolo-tlon  function  directly  as  an  output  volto^i 
which  coiild  b©  moaourod  and  recorded  by  properly  oolectod 
Inctnaaents . 

Th©  blodt  dlagroia  of  such  an  anaxoG  computer  do  signed 
and  develoixsd  by  Kennotli  V/.  Ooff"^  is  alio-wn  in  Fig.  2.1. 
A  photoG^aph  of  the  correlator  le  shoim  In  Fig,  2.2.     For 
^ba  puri:)0oes  of  r.naly.^dG,  the  eoaimtor  lo  ccnpocod  of 
tliroo  separate  units 

1,  tlae  delay  nygtem 

2.  multiplier 
5«  Integrator 

li&Mx  of  thoce  unitcj  raay  be  operated  indeixjxidently  of  the 
reaainine  conpomonts . 

2.2   (a)      THE  HTME  DiCEAY 

The  tlae  tielay  moohanlsm  employed  in  the  caiiputor 
being  considered  conriots  of  aii  electro-mechanical  dovlc©, 
^S»  2.3.      -\s  Bho\«i  in  the  block  dia^ran,  the  t*»ro  slGnal 
voltai:;eD  boinc  uoed  B.Te  fed  Into  nej^rate  channele  each 
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Fls.  2,2    mdtogmi&i  of  aimlos  ooiroXatlcn  comr^ttter 


If 


TRACK  1  (FIXED  DELAY) 


ERASE 


INPUT 


RECORD  AMP 
AND  PRE 
EMPHASIS 


-»-" 


RECORD 


RECORD 


—      ,  ^~N'^(P0SITI0N 

-t-T'.n^'^  /for  max 


DELAY) 


OSC   AND 
POWER  AMP 
FOR  BIAS 
AND   ERASE 


CHANNEL  1 

REPRODUCE  AMP 
AND  POST  EMPHASIS 

CHANNEL  2 


RASE 


f(t+-S) 


OUTPiJT 

► 

f(tt4) 


REPRODUCE 


r,w'^    (POSITION  FOR 
—  ^        _-     ''      ^  ^MIN    DELAY) 

\  j;'^   RECORD 
TRACK  2  (VARIABLE   DELAY) 


ffUh  d«9  n^m  t$mfm  ^  memUi$  wm  tmu^ 


18 

having  Its  ovm  variable  Gal«  ampllflGr,  raacnotlo  record- 
log  head,  mtignotlc  rocorcllnc  track,  playback  head,  erase 
head,   and  playback  ampllflor.     Th©  record  head  of  channel 
one  lo  fixed,   viioreas  tho  rccor(5,  head  of  bhannol  tv/o  Ic 
movable.     If  tho  trancit  times  frcsn  tho  record  to  the 
playback  heads  for  a  point  on  tho   ourface   of  tho  ro cording 
drum  aro  equal  for  both  cShannelo,   then  no  rel-^/tlvo  time 
delay  \dll  be  introduced  between  the  two  Input  voltaseo. 
If  tho  transit  tiae  between  record  and  playback  heads  is 
greater  on  chanr^l  two  than  on  Giuutntl  one,   then  a  rela- 
tive time  delay  equal  to  the  difference  in  trnnalt  tlnws 
is  Introduced.     For  Identical  Input  voltages  to  each 
ifiMmnel,   tho  oxxtput  of  diannel  cme  would  bo  f  (t)   and  the 
output  of  channel  two  vrould  be  f  (t-T')   vjhore   T    la  tJie 
difference  in  transit  time* 

Thlo  relative  time  delay  between  channels  aos  and 
two  can  be  varied  contlnuouoly  from  -15  mlllinecondQ  to 
190  mllllr.ecc2ids  by  laeano  of  s^ar  train  and  dialn  drive 
synt^a  poRfored  by  a  low  speed  synchronous  motor • 

The  frequency  responoe  of  the  oyotem  Is  flat  to 
ifithln  -  2  db  over  a  freq-uency  ranso  from  ICC  cps  to 
10  kops, 

2.2   (b)      TfIS  !4UL!rinLIi:R 

The  multiplier  companont  of  tho  correlator  is  a 
quarter-dlfforeaice-Gquarlns  device  x-fhich  reduces  the 
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prooens  of  multiplication  to  that  repr^r^nted  by  the 
relation 

fx  f«  =  1   (fx  +  f^)^  -   (fx  -  U)* 


i 


The  process  of  producing  a  voltage  proportional  to 
the  square  of   (fx  +  fa)   and   (fi  -  fa)   Is  per- 
formed by  a  sqiiarlng  circuit  connletlnc  of  two  6B8 
pentodos  operated  v;lth  grids  In  push-pull.  For  opera- 
tion In  the  region  where  the  trannoonductance  and  grid 
voltage  are  lineally  ixjlated,  the  plate  cixrrent  for  one 
tube  can  be  approximated  by 

1-  =  a  -f  b©^  +  ce„*  ♦  de-^ 
o  6     5      S 

For  the  two  tubes  operating  in  parallel  vdth  a  common 
plate  load  as  shown  In  Fig,   2.4,   the  current  throu^i  the 
ccraraon  load  Is  the  sum  of  the  Dlate  currents 


1.   =:  1^4    +  1^3  =5  2a  -f  See** 

L         pi         p  1 


lihe  output  voltage 


p 


^o  =  Sb  -  ^L  ^  =  ^  -^  ®*l" 


If  two  sq\aaring  circuits  are  operated  in  push-pull 
and    e.     is  made  equal  to       (fa.  +  f«)       In  one  sqiiaring 

and  equal  to       (fx  -  ^a)       in  the  other  squaring  clrciilt, 


00 


JPjkj^    QJt^::iphf\ed   -.■hematic  of  the   ^.i\  larinf^  <  ircuit. 
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then  the  outpfut  voltage  will  be  proportional  to  the 
product     fi  fa 

Variability  In  the  pe rf orsnrince  of  the  correlator 
generally  results  fi»om  vmbalanced  oondltiono  In  the 
multiplier  clroulta.      'hen  the  miiltlpller  la   oporatlag 
oorroctly,   the  output  voltage  should  be  zero  vAien  there 
la  no  Input  and  vjhen  there  l3   only  one  Input.     "JSiono 
balanced  conditions  correspond  to  the  requlremontG  that 
QjcO  =  0  and  0x1  =r  0  =  1x0. 

The  balance  for  0x0  =  0  la  made  by  adjuctlng  tlio  d© 
bias  of  the  four  6b8  toabeo  for  eqtml  quiescent  pinto  cur- 
rents v/hereas  the  balance  fbr  OkI  =s  1x0  «  0  is  made  by 
adJUGtlng  the  ac  levels  applied  to  ea<ai  6D8  tube» 

Much  attention  waa  focused  upon  the  design  require - 
BMRte  necessary  to  reduce  do  drift  to  a  dOGroe   that  would 
permit  continuous  operation  for  periols  up  to  two  to 
three  hours  without  adjuntment   • 

2.2    (c)      THE  INTEXJRATOR 

Integration  in  the  cosaputer  is  accc»!iplir>hod  by  means 
of  either  a  low  pass  filtering  system  or  a  stepping 
IntGsrator.     For  acoustical  laeamireiaentG,   a  simple  RC 
low  mas  filtering  clrctAit  Is  satisfactory.     Hhe  corre- 
lator has  oeveral  su<di  circuits  vrith  RC  tliaes  of  C.5, 
1#0,   2.0,  4,   8,   and  16  neoonde.     Selection  of  the  proper 
RC  time  is  another  important  factor  sovomlng;  optim\£a  use 
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of  tho  cruli^raent. 

Tlie  6e  cmtmit  of  th©  RO  clrciAlt  Is  amplified  find 
then  r>i!c»eed  aororp  the  In^ut  of  a  dc-ac  choD-oor  v;h.lch 
convert"   tho  '^c  voltfic©  Into  a  pronortl^viai  AtHj  c^b 
signal.     Tlilo  AOO  ods   signal  Ifl  tlien  am'^llfled  -^pialn  and 
r>lp.ood  acrons  the  ovitput  termin£?.lo  of  the  intecrator. 

©le  400  ops  algnftl  nay  ^  neapiired  directly  :flth  a 
vaouum  tube  voltraeter  or  an  oncllloGCOt^e,   or  It  may  be 
rectified  and  rocordod  by  meane  of  a  do  graphic  aometer. 
A  llnecir  plot  nuch  as  is   obtained  on  a  f5raphlc  aanaeter  la 
extreaely  UBOf\il  Rinoe  it  reprenents  tho   shape  of  th© 
COTiputod  correlation  function.     Information  regarding  the 
relative  amplitude  in  decifoela  of  suoceoQlv©  peaks  of  tlie 
congelation  curve  la  desired.     3uch  information  may  be 
obt^ained  by  moans  of  a  logaritiunlc  graphic  level  recorder 
^Ich  records  the  correlation  curvos  directly  in  db. 
When  used  together,   the  liiteor  and  loGarltiimlc  plots  gIv^ 
a  continuous,  detailed,  visual  record  of  the  shape  and 
magnitude  of  tlie  correlation  curves  as  a  function  of  time 
delay. 

All  of  iiw  information  preaented  in  tiie  foregoing 
dlsoxiSGion  of  the  analog  oorrelaticaa  computer  has  boen 
talcen  from  the  doctor s»>l  thesis  of  Kenneth  ./.  Goff  • 
Only  those  points  vAiich  are  ccinoidered  neoessary  for  a 
basic  understanding  of  the  computing  proco^s  are  pre- 
ewited  here.     For  a  mare  detailed  treatment  of  the 
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ddslGxi  roqulremont  and  the  porforaano©  charactorlntlcn 
of  the  coraputor,   tho  reader  Is  roferrod  to  either  rofor- 
onoe  3  or  to  the  recontly  complGted,   unpubllohad  maniial 
"Instruction  Manual  for  Analog  Correlator"  by  Kenneth  ■<. 
Goff,   -luguct  1954,  H.I.T,   AcouotlCQ  Laboratory. 

2.3     DATA  GATHERING  ACCE330RISS 

ITie  aooossorles  selected  for  G^tlioifSK  "^^"^^  ^^^  ^-^ 
f ollowB : 

1,  The  nolso  source  consisting  of  an  ©loctronic 
noise  voltage  Generator,  a  pov;er  ampliflor,  and 
an  Altoc  728b  loudspeaker  in  a  clocod  box  baffle, 

2.  The  mlcrophcaies  used  were  Altec  Lancing  type 
21  DR-150  vdth  157 -A  bases  and  P518A  Altec 
Lansing  pov/er  supplies.     The  response  cuirves  for 
these  mlcroi^ones  are  flat  to  within  «  2  db  for 
a  frequonoy  range  of  froa  100  cps  to  10  kcps. 
The  voltage  amplifiers  ur.ed  with  these  micro- 
phones vfore  wide  band  amplifiers. 

3*   Since  the  reducing  equipment  ^me  not  portable, 
provision  for  rooordlns  data  was  required,     Aji 
Ampex  model  350-2  twin  track  magnetic  tape 
recorder  (Fig,  2,5)  vras  uoed  for  this  pixppoce. 
The  particular  apparatus  used  had  a  frequency 
respanr.o  flat  to  v/lthln  -  2  db  over  a  ran^e  of 


s» 


ns*  S»5    aiotoGraiAi  of  mmx  nagnjitle  tape  r»«oinl$r 
«lth  a  tape  loop. 
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150  cpa  to  15  kops.     Tills  oqijdpiaDnt  vras 
portable  ajid  yot  capable  of  roprodxiclag  a 
r«lAtiTely  dlctortlonleaa  replica  of  tho  In- 
put voltage* 

2,4     ACCE:>30RIE3  for  PKS3EHTma  RESULTS 

The  corrolp.tion  function  an  compnbed  by  the  arj?.lor> 
correlator  1b  presented  In  tlio  form  of  r.  4oo  cpa  voltas® 
and  a  roctified  dc  voltage.     It  Is  therofore  possible  to 
tt»as\xro  and  record  the  correl^vtlon  fimotlon  by  aeane  of  a 
variety  of  electrical  meaoiirlnG  InntrrciontB. 

lUxe  aooeasorlos  used  for  xa»aBtirlnc;  volttices  are: 

1,  Model  3CX}  BallantiHG  vacuiaa  tube  voltmeter. 

2,  Typo  208  Du'-iotit  Catliodo  Ray  OGcilloacope, 

Perraanent  i^oconls  of  the  orutpitt  voltage  w»re  aad» 
by  meanB  of: 

!•  Modol  AW  linear  scale  3:3terllne  Anf^us  Grephlo 

MlXllaxmieter. 
2,  Bruel  and  Kjaer  log  Gcalo  Grar^lc  T^vel  Recomer. 
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III.        TRANSIEWT  RESPONSE  OP  ROCMS 

3.1      INTRODUCTION 

Many  studies  havo  been  made  on  the  oubject  of  oound 
flelda  In  rocmB.     Genorally  there  stiJdlGs  have  been  based 
on  one  of  two  distinct  approaches  to  the  problem  of  room 
acoustics : 

a)  a  study  of  the  effect  of  rootriB*  boundaries  on  tlie 
sound  field  aacoclatod  with  steady  state  pure  tone 
aignalfs. 

b)  a  study  of  the  effects  of  rooois*  botmdarles  on  the 
QovnQ&  flold  associated  v/ith  transient  signals. 

Therefore,   in  attoraptlns  to  obtain  a  ccmplete  qtumtl- 
tatlve  ansessnient  of  the  effects  of  boundaries  on  a  sound 
field,   one  mi£^t: 

a)  excite  tho  room  with  a  non-dire ctlcmal  steady  state 
sinusoidal  slgjial  and  map  the  amplitude  and  pha«« 
relations  of  the  sound  preoaure  for 

1.  all  freouencies 

2.  all  possible  observer  positions  in  the  room 
3*  all  possible  source  positions  in  the  room 

b)  excite  ihe  room  tflth  a  pulse  of  extremely  short 
dtaratlon  and  measure  the  rel^.tive  amplitude  and 
direction  of  propagation  of  th©  associated  vrare 
fronts  for 
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1»  all  time 

2,  all  obserTor  positions 

3«  all  oourc©  positions , 

If  on©  has  all  of  the  data  aasoclntod  with  either 
one  of  the  above  mentioned  experimental  apnroicheo,  ho 
e«a  mathematically  deduce  all  of  the  InToimritian  for  th« 
other  method.     For  •XM^e,   on©  can  synthesize  th©  ptilo© 
used  in  terms  of  a  Fourier  series  of  sinusoiciale,     Th«n 
for  each  frequency  component  one  can  pick  the  relative 
ampllt\ide  and  phaae  angle  of  the  sound  preosuro  at  any 
point  in  the  room.     By  \fBi5l1tinG  thon©  amplitude  and  phB.Be 
itteasurements  in  accordance  ^flth  the  Fourier  analysis  of 
th©  pulse  and  then  adding  the  effects  for  each  frequency, 
the  complete  responre  of  the  sound  field  to  a  tranoiont 
puis©  can  be  c^educdd, 

on  th©  other  hand  on©  can  visualiz©  a  model  of  a 
room  excited  by  a  pulne.     In  ouch  a  model,  the  walls  of 
th©  rocjm  would  b©  eliminated  and  their  effect  on  the  sound 
field  would  be  simulated  by  an  infinite  multidimonolonal 
array  of  image  source c,  eac5i  of  \4hich  is  located  in  a  cell 
of  th©  sam©  sis©  and  shap©  as  th©  actual  rocnn.     Eadh  imag© 
source  is  positioned  so  that  the  transit  time  of  th©  sig- 
nal from  the  im^-g©  source  to  the  actual  obBorvor's 
position  entials  the  transit  time  for  the  signal  arriving 
at  th©  observer's  position  by  a  particular  sequence  of 
reflections  from  the  room's  botaadarlea. 


I 
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v/lth  thlo  ma^-el,  one  can  now  Ajjeclfy  that  oach  non- 
dlrGCtlonol  source  Is  In  phas©  and  is  radiating  a  froe 
travellnc  olnuDoidal  aignal.  By  adding  up  the  contribu- 
tion from  eacfii  image  source,  one  con  nap  ths  amplitude 
and  phane  relations  of  a  room  excited  by  a  steady  state 
source. 

The  above  dlocussicn  assumee  that  the  boundaries  of 
the  room  are  Infinitely  hard  so  that  no  phane  shift  or 
decrease  in  amplitude  occurs  vjhen  a  signal  is  reflected 
from  the  vmllo.  This  ass\imptlon  simplifies  the  problem, 
but  tlio  reoultlns  calculations  do  not  conform  v^lth 
reality.  It  is  necessary  therefor©  to  connider  tliese 
phase  Ghlftn  and  attenuations  occurring  at  the  boundaries 
of  the  sound  field.  In  the  transient  treatnent  of  sound 
fields  one  may  account  for  the  boundary  effects  by  as- 
sisning  to  each  bounding  surface  an  Impulse  roBponne 
H(a)) ,  In  studying  the  steady  state  behavior  of  bounded 
soTJnd  fields  an  absorption  coefficient  and  a  phase  shift 
are  aasisned  to  oach  v/all. 

Any  such  complete  analysis  of  a  bounded  sound  field 
would  be  extremely  tedious.  A  derire  to  be  practical 
v/ould  incline  one  to  restrict  tJio  analysis  to  a  degree 
«W9BBienGurate  vjlth  the  amount  of  detail  nececaary  for  the 
proper  utillaation  of  the  room.  In  most  rooms  we  are  in- 
terested In  producing  an  acoustical  environment  conducive 
to  s^oi  atKllo  communications.  Therefore  \m   may  Generally 


29 

ro strict  our  connlddrntlfm  to  the  atidlhlo  froruoncy 
raai$».     Audio  coniraunlc'^tlans  v/lll  sonorally  fall  Into 
ome  of  threo  cr-toGorlcs. 

1.  two  way  vocal  coraTmmlcntlons  b©V.^en  Indlvldualt 
?•  Toeal  eoBBBoni  oat  Ions  bot>^an  <xo&  epeator  and  an 

audience 
3»  inotnMental  coeamuni cations  "bet^^on  muslciojis  and 
an  audience 

Traditicanally,  each  tyi:>e  of  comriiunicatiGin  montioited 
above  has  associated  with  it  a  cei^ain  phyalcal  arransement 
of  Goiirce  and  obsorver.     If  \^  accept  the 00  traditional 
arrangenonts  ac  being  charaotorlQtic  of  the  type  of  ccra- 
mimioatlon  to  be  unod  in  tlie  room  then  ^«  may  further 
limit  our  analysis  of  a  bounded  sound  fiold  to  th08« 
source  rjid  observer  positions  of  practical  interest. 

For  practical  purpoeofl  wi  may  further  restrict  otrr 
analysis  by  cons  id©  ring  only  thone  types  of  scrundc  cliar- 
acteristic  of  musical  or  vocal  sources,     Cudi  simplifioa- 
tionc  furtlaer  i^strlct  both  the  frequency  coin":>onitiQn  of 
the  toe  tine  signals  and  pulse  diiration  and  ptalse  inter- 
vals for  tranclont  signal  tec ting. 

Up  to  this  point  \^  hiive  boon  din  cussing  the  tv/o 
general  ai^proaches  to  the  study  of  room  acoustics  and  the 
various  simplifications  that  may  be  employed  in  this 
study.     Those  two  approaches  have  been  quite  thorot;i^ly 
explored  and  the  restating  body  of  data  coupled  with 
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subjoctlve  evaluatloIUl  of  tlio  acoustical  envirozjmont  In 
roomo  hao  led  to  th©  conelunion  that  varlationo  In  both 
the  Dteody  nt-^.to  charp-ctcrlftlcs  and  the  tranoiont  diar- 
actorlstlcs  of  oound  flolde  account  for  dlfferoneoa  In 
thi  acouGtlcal  quality  of  a  room, 

W,C,   cablne     re  colonized  the  Import  nnce  of  the 
transient  behavior  of  aound  fields  and  attenpted  to  re- 
late It  to  a  subjective  evaluation  of  the  quality  of 
rooms.     As  a  measure  of  transient  reapons*,   Sabine 
defined  a  '^revorberatlon  time"  ^ig  lenpith  gf  time  for 
Hit  mim  square  pyost^ur?  St  a  e\4tab3,y  ghoaeq  d4.ntrlb^- 
%Xm  2t  ^PW^  B^v§s  tQ  6X^XT^%^h  lo  ^m.  ?ii;3,;34onth  ^  XM 
orlr^ln^il  Intensity,     Re'ferboratian  time  as  defined  by 
Sabine  is  a  moar^uro  of  tho  rate  of  do  cay  of  the  aound 
preanureo  caused  by  a  steady  state  alnusoldal  signal 
lihich  la  auddenly  discontinued.     Various  criteria  "based 
on  reverboration  tlso  vs  tho  slae  of  a  room  ^^Idi  la  to 
be  used  for  a  apeciflc  application,   and  revorboration 
tJUaft  va  frequency  hr^ve  be<5n  devised'  ,     At  present, 
design  ctirvca  of  this  type  offer  the  most  satisfactory 
approach  to  the  construction  of  rooma  having:  good 
acoustics  and  the  treatment  of  rooms  having  poor  acous- 
tical characteristics.     Hov/ever,   it  Is  recognized  that 
acoustical  denlgn  critorla  other  than  reverberation  time 
SKbBt  be  considered  in  the  correot  design  of  rooma*     For 
instance,   one  research  rtroup     s\iGS®"ted  that  the  flra% 
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20  db  of  aeoay  of  a  sound  In  a  rooea  lo  of  prlaary  Im- 
portanoo  In  <31fferentlr,tlnc;  botvroon  two  rocaaa  vrtiloh  hav© 
approxlmr.tely  equal  ovrjr-all  rovorborcitlon  times, 

n.A,  Mason  and  J,  Ilolr^,   in  the  courro  of  t2\Glr 
otudloB  of  the  acoustics  of  theatoro  by  ■•ans  of  puln©g 
or  "toffie  buroto,'^  Intimated  that  the  time  and  anplltude 
dlotrlbutlon  of  reflected  tono  btirotc  mlsht  servo  as  an 
aid  In  ©valuatlns  the  acouotlcal  quality  of  rocrao. 

In  the  Introduction  to  "Pulse  L^tatlotlco   'ooalyslQ  of 
Baoni  Acoustics"'  the  authoro  state,   "...  Vheoe  facts  Incll* 
cat©  that   ...  tiie  flmt  fev;  reflections  are  primarily 
responolble  for  cortnln  Important  features  of  tlie  acous- 
tical character  of  rooms   ...     It  should  therefore  be 
worth\^?hlle  to  ntudy  this    'short  texm*  tranolent  respofio© 
by  a  method  of  Images  In  which  all  of  tlie  wave  proportlos 
of  the  Imaoe  so-urces  con  bo  considered   {l,e»  vrhore  the 
assumption  of  an  incoherent  source  is  not  made).     Further, 
if  an  imago     array  satisfying  the  bourdary  conditions  can 
b©  fousxi,   one  should  be  able  to  treat  this  array  statis- 
tically and  thus  obtain  the  long  terra  average  transient 
responco  as  \mll  ..." 

As  indicated  by  tJie  above  quotation,   seme  method  of 
•©asurins  tho  transient  rosponn©  of  rooms  is  required  for 
an  evaluation  of  the  effects  of  boundaries  on  a  sound  field. 
However,   ccoe  estlmr.te  of  tho   nteady  state  behavior  of 
roans  is  equally  desirable.     Cne  exi>erlmental  raothod  of 
aeasurins  the  transient  response  of  rooms  has  been 
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ougs^oted  whloh  Boerao  to  canbiLti©  Bom  ef  th©  i>ropertloa 
of  both  otoody  Bt-^te  methoae  and  pulce  motliodE}-  .     This 
aothod  employa  cori^l^itlon  tooimlquca. 


3.2    TIC  cor.  si/-.ti-:n  '^;i\Lx:;ic;  07  T?JVKJi;i:KT  iggpcsios 

Now  that  we  have  su^cooted  oooe  potentially  useful 
mcthodo  for  maJdUig  acountlcal  measiireraento ,  It  would  be 
wise  to  exploro  briefly  tlie  mathooatlca  of  thlo  method. 

It  can  be  shown  dlx^ctly  that  tho  corroct  applica- 
tion of  correlation  moa.Burlnf$  toclinlques  to  bounded  sound 
fields  -vili  yield  tiie  Impulse  response  of  that  bounded 
Bound  flold^, 

neforrlnc  to  Pl6»  ?»1,  ym  assume  that  tlie  nols« 
Dourco  I3  IcllvDirlnG  a  slsiial  having  the  coni:^le:c  spoctrum 
P_(a))   to  the  tvro  channels  of  the  correlator.     If  tlie  net- 
works  IxaVG  transfer  functions  Hi  (3)  and  H«(3)   reGi:)octl7ely, 
then  the  outpiut  of  the  notv/orks  are: 

Pi(o))   =  Hi(jca)   ?g(«)  (3.2-1) 

Ps(a3)   =  H.'.CJcd)   Fg(«))  (3.2-2) 

The  cross  poi.'er  spectrum  for  Fi(<o)  and  r«(a))  is 
*21^®^  =  l^^  2  [Hi(ja)  Fg(o)  Ha(-»  F^)] 

=  ^ne  (ca)   Ha(ja))  Ha  (-jcd)    (3.2-3) 


■ 
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5* 
where 

Is  the  poorer  denoity  npootrtra  of  tlie  nolo©  scniTM* 

Hovr  the  eorrelaticn  Is  the  Fourier  MalMi  tMMKl 
of  the  poft^r  <5enfllty  spectrum 

(p^,    (r)  =  |#„„(c»)  Hi  (Jo)  H«(.jo)  e^^^dm 


(3.2-5) 

From  Fig,  3,1  \m   can  nee  that  frequency  variation  In 
the  source  can  be  considered  as  common  terrae  to  both  ohan- 
nalc  ono  and  two.  This  factor  then  con  be  conoldered  aa 
a  eonotant^  l«e* 

<^  21  (  ^)  =    $  as   )  ^^^^^®>  Ha(-jG))  e  ^^^d» 

(3.2-6) 
If  mi  nmt   conDldor  th»  epecial  case  ^4iere  H^  =  1 
and  y©  \flsh  to  auto-correlate  the  output  and  input  of  th,e 
system  having  a  transfer  function  of  lUCjca),  the  cross- 
oorrelaticn  function  la 

^  21  ^  ^>  =   ^  8«    f Ri(Jcd)  0  ^^^  dta  = 

2»T   ^  gg   hi  (t)  (3.2-7) 

^hus,  the  cros6-oorr«'latlon  function  Ib  exactly 
proportloaml  to  the  impulse  i^oponoe  of  the  syotem  having 
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th©  trancfor  fimctlon  Hi(jo).  If  the  channel  harlaB  the 
transfer  function  Hi(J(o)  lo  a  roQci»  then  the  correl:\tion 
function  Is  proportlcmol  to  the  Impulse  roeponse  of  the 
rooo* 

H'o  hav«  tallied  a  little  about  the  nature  of  tho  ocjfr- 
ftlatlfwi  function;  let  un  novr  conolder  Its  flliApe  for  a 
correlr.tlcn  analysis  of  the  sound  field  of  a  room*  For 
this  purpooo^  \m   shall  firot  define  a  fow  tome* 

Usinc  Laplaoe  Tranofonnation  we  made  the  definition 


Fx  (cd)  i  1    \    fi  (t)  e  "^^  dt     (3.2-8) 
2n 


whence 

tx   it)     a      I  ^x   U)  e  '^'^^   do 
fbe  Multiplicrtion  Ihoorem  states: 


I 


«0 

jot 


air  J_^^i(to)      Fe((ia)  e  ^^    do     =s 

(t)  fa  (t-^)  dt  (3*2-9) 


Now  dividinc  ©<*^ch  cid©  by      2T      'tfdsezv    T    =»  tiiae, 
axxd  then  taking  the  liadt  as     T-*«* 


j     Fi(oj)  P, 


lis  StT  I  *    ^ 

T">-»  2?  I     Fi(o)  IvTo)  e  ''^^  do  =5 


T^*  If    J-T       fi(t)   tf>A%-r)  dt  ^(Pr,^   (T) 


(3.2-10) 


In  -Uio  application  of  correlation  todmlqnea  to  tNi 
troaslotnt  rospfina*  of  roonm  fx(t)  la  tho  roltoQfi  produced 
at  the  loudcptalBir  t^xnlnale  by  a  vld«  load  noloe-voltaoi 
tfmBratOT  and  f«(t)  Is  the  tt>lta^  output  of  a  mlnrophcni 
placed  at  aoa»  <^eervutlcti  pl&ot  In  thA  sound  field. 

The  voltaet  produced  by  the  microTiiono  is  foCt),     It 
conolntn  of  the  sum  of  fi(t)  nodiflod  by  tli©  tmncfor  (diar- 
aoterictic  n  of  the  Iqudap^akor  and  mLcrovtioane ,  plun  a 
oiwsaatlcai  of  rofloctiotts  each  of  vihich  oquals  fx(t) 
delayed  by  a  tino  equal  to  the  path  Ittngth  of  trojicDiQeloKii 
If  divided  by  the  opeod  of  sound  a;  and  modified  by  H  and 
A  tawmr.for  characteriRtlo  B  associated  with  its  tranamio- 
sicai  patli. 

The  fro  iusncy  daaaln  ropre Dentation  of  euch  c,  -/oltage 
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jo  iVo  jo  li/c 

F«  s:  PjL  H  B^e  "^     +  Ft  n  Bi  e  + 

jo  In/C 
PiHBa  e  + (3»2-ll) 


B'o  are  complex  ooefflclento  whidi  account  for  the 
iJiaee  shlftc  at  the  reflecting  aurfaoesi  the  attwuiatlon  in 
aaplltude  acaoclated  with  Inverse  squar©  dlnruniticsa  of 
Wdiatod  sound  IntefXieity^  and  the  attenuatlcai  in  *arr3lltude 
occur r ins  at  the  rofloetlnc  surface. 


B    ^    JBJ       e       ^^^ 


<^  -Jco  I/O 

t^«  »       ^    ti'n  H  O  ^         3^  (3.2-12) 

s  n  pamrr  Awadty  apaotrua  =    lla         tt  lF(o))i  J 

(3.2-13) 
theai  It  follc^rc  that 

f'(«»)  =ito  /t.s        o     ^^  (3.2-14) 

T^  '    tr 

P(a))  a         /f(o)/  •     1^ 

!fow  stibtstltutins  la  th©  multiplication  thoawwaj 

C^         g  J   (eat  ♦*=,  -*x)  ^       (3.2-15) 


c 


>?hoii  vjD  corrolate  P*  and  ?«  ualns  ©qur.tlon  3.2-15 
and  3.2-12 


.  ^  (Pn-'^V")] 


dcD 


3« 

(3.2-16) 

Sqiiatlcsi  3  •  2-16  then  Giv«8  the  crose-con^latlon 
<P*2(7-)     betwMm  th»  Input  rolttx^e  to  tliQ  lowdBpeolatr 
and  th©  output  volta@s  of  the  microphcqno  In  toiwc  of  th« 
pmf&r  dfanalty  epectrun  (s)   of  a  nolDe  voltoc^t  emmrator, 
tlio  oomblnod  tranofor  function  (H)   of  tho  loudopoalcor 
aal  nlcrophcaiD ,  and  the  coraplox  tran«f\ir  function  of  th© 
tronnnl colon  path,     This  latter  trancfer  functlott 
/y  o  3^      a«ocrtb..  th«  .cmplax  »ollflcaUon  In  th. 

fiomsfl  pirosQui^  aunplltudt  eptetrm  caused  by  Itc  traveling 
th0  distance  1^^  Involvlns  ^ny  number  of  reflections. 

For  the  idealised  ease  M&wvet   (a)  souroe  radiates 
nolee  of  a  flat  spoctrua  frcsa  ax  to  <»o  and  zero  outolde 
of  this  apoctnm;   (b)  the  transfer  responoos  of  the  loud- 
WptaSsoT,  mlorophosnea^  and  walls  qX9    XIQ"  $  It  has  been 

TO 

Shown      tJiat  tlie  croos-correlatlon  ounm  for  tho  loud- 
8pea;tor  and  siicrophone  voltasos  14.11  have  equal  peAloi 
for  Ga&i  point  i^ere  the  tlae  delay  oorrospondo  to  the 
tranclt  tlsie  fi'oia  the  aouroo  to  the  microphone  and  that 
the  cross-corr»olatlon  curve  ^dll  danp  out  to  a  eabSU 
iralup  beti^oon  those  points.     The  rapidity  ifith  which  It 
dsfips  out  'vfill  depend  upon  the  bandwidth  (oa  -  oi)  and 
the  term  (o«  4-  ox)* 


fi 
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With  tha  f«mgfliiig  aAtlMHBo.tlcs  In  mini  w»  aajr  nov 
visualize  ::  room  model  that  will  he  uo^Tvl  for  corrala- 
ticox  analyalB.     This  nod*!  will  \m  exactly  the  same  ae 
liie  Modtl  Mnii«ie4  prerlouoly  in  seotlon  3«X  In  connec- 
tion v/ith  iiulae  analysis  of  rootts.     In  this  aodtl  tlie  walla 
are  roplaood  by  an  array  of  laaee  eotircoo  ea^  of  Which  ie 
located  at  a  dietanee  equal  to  the  lensth  of  the  trettttBlQ- 
sion  path  of  the  rtfleeted  sound.     The  outinit  of  eacii 
souroo  day  be  modified  in  aocordanc©  with  the  nuabor  of 
Inflections  thr-t  ooetir  in  the  path  of  transmlsaion  from 
Hie  uwtrm  to  the  db«€^t<i^r*e  position. 

In  thoir  dlGousslon  of  pulae  static  tics  theory, 
B«lt»  Doak,  and  •festervolt    dovolopod  the  matheafttlCtt  of 


a  model.     In  their  model  Miflh  Isa^i  source  occupies 
a  separate  iasgo  replica  of  the  physical  rocaa.     tlmm 
iMMSe  rooms  or  cells  arc  designated  by  three  !iuii^bev« 
{It  a,  n)  eacSi  of  which  can  tato  all  lntejTa.1  valtjes  from 
rainuc  infinity  to  plus  infinity.     I^o  notation  deoisnating 
the  tm&^  sources  is  also  an  Indicator  of  the  number  of 
refloctions  that  the  slsnal  undorgoos  in  traveling  from 
tl^  GOurCG  to  the  obsorror, 

^  iShall  now  aal»  a  few  simple  adaptations  of  this 
spatial  model  to  aid  us  In  visuallains  tlie  corrolation 
analyclo  of  the  problem.     All  of  the  sourtjos  oliall  be 


AO 

grouped  in  aooordonot  t/lth  the  nunber  of  rofleotlons 

that  occur  in  transit  from  tho  acfurco  to  the  obeerTor. 
TtiUB  tlio  oooireoB  uhooe  sound  io  refloctod  tvdoo  In  trans- 
it are  doslcsnated  as  ooeond  ozxlor  ima^  souroos. 

fThe  nuAber  of  SjmmP  eoixroeo  in  ea6h  ordsr  depends 
tipcn  tho  rocaa  (sooaotry.     For  oxanple,  tlie  aodol  for  a 
ivetone^^^  Tocm  oan  be  ^ctom  doun  Into  Its  loo^  com- 
ponmte  w  foiioifst 


fr     ■  ■ 

fiSSB 

Actual 

1 

Ist  Order  2^6* 

6 

2nd  Order  Zaaoi 

xa 

3rd  Order  laaae 

3S 

etc. 

An  outline  of  tho  valine  onolortKI  Ijjf  tte  actual  plus 
tike  six  first  order  iaa^a  souoroes  for  a  roctanBular  rooa 
is  shoim  In  ?ls*  3. 2*     The  number  of  liaae*  sources  of  ea^ 
order  can  bo  oooiputed  froo  the  formula 

Hisabor  ist  2  ^  ^  -¥  B       21.      k  (3*3-l) 

iidxere    K    Id  the  order  of  the  reflection  and  k  =  (N-1)  • 
For  example,  for  the  fourth  order  array,  there  will  be 
2+16+8  (3+2+1)  »  66  Image  sources. 

Other  factors  of  laportanc©  in  correlaticn  analysis 
are  the  distance  beti^oen  Ins^  seyiircos,  and  t2ie  distance 
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Plg»  3.2      Volvcao  enclosed  by  actual  and  first  order 
iTiage  sources  for  a  reotansuln.r  roomp 


i 


from  tho  laa^i  souroea  to  the  obsorrar's  position. 
TbOQ%   faotors  may  bo  oanputed  from  formulas  9  «nd  11 
of  roforonoe  6* 

nit  «dY»ntas««  gainod  by  roferanoe  to  this  spatial 
aodol  of  a  roo9A  are  threefold; 

a)  tho  system  givos  a  visual  prosentaticn  of  tbm 
problsms  encountered  in  spaoe  or  tlms  flltmwiMfjt* 

b)  the  results  of  the  correlation  onalyaic  may  be 
compared  '.flth  those  of  other  reaoarchers  who  used 
pUlae  nethods  of  analysis. 

c)  the  mathom-ntios  for  puloe  analysis  developed  by 
the  MBB   of  this  model  may  be  applied  directly  to 
oorrelatiasi  analysis. 

The  expertnewtal  attasuriments  of  the  transient  re- 
sponse of  rooms  were  taken  by  cross  correlating  between 
the  input  voltaLge  to  a  lofiidspeaker  ax»l  the  output  of  a 
microphone  located  in  the  room.  A  hard  plaster  room. 
Room  10-390-A,  ifoich  had  been  previously  analyzed  by 
poise  skethods  vas  selected  for  the  preliminary  terts. 
Both  the  source  and  microT^one,  as  vjoll  as  tho  teat 
positions  in  the  room,  vnsre  solocted  to  give  results 
t«a;»arable  to  those  obtainsd  by  Bolt«  Doak,  and 
v^stervelt®. 
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Source  Position  J 


The  r#0tQ.t8  for  on»  p&rtloxolar  test  rm  aro  shoMn 
In  ^is*  3  •3*     shoim  also  aro  the  ecntparable  rastslts  of 
a  pulse  analyolo.     In  this  run,   tho  Boairo6f  a  Model  173 
Janes  Tianelng  horn  driving  unit  vdth  a  Model  1217*1290 
Janes  Lanalng  throat  and  leno  assembly,  was  plaoed  in 
the  comor  of  tho  rooan  facing  the  wall.     ITho  mlcro^^hone, 
an  Altee  21*BR-150A  was  looatsd  13  foot  from  the  ocurce 
alone  t^e  ohort  wall  of  tho  room. 
The  portinent  dlmenBiane  aro: 

Rocd  diisensions:       length      Uc  :=  23  IViot 

width        Ly  =  13 •^  feet 
heii^t      L2  =  8.4  foot 
X  =  1.0  inohes 
y  =  1.5  inchos 
z  =  C.5  induMi 
u  =  1.0  inchos 
▼  =  13  toot 
V  2=  1.0  inches 

Th9  teot  signal  was  rendom  nolee  generated  by  an 
eloctronlo  noise  voltage  genontor  and  aaB5llfled  by  a 
XOdel  20-W2  HolntoBh  amplifier*     k  1/3  octavo  band 
analysis  of  the  loudspeaker  volta^  and  the  raiorophono 
voltage  is  shown  in  Fig,  3.4.     Bie  loildspealosr  volta^ 
and  the  aicrophono  volta^  wore  recorded  simultaneously 
on  separate  trac&s  of  a  twin  track  tape  recorder.     In 
reducing  the  data,  a  eecticai  of  th©  recording  tape  was 
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out  cmxl  m»a»  Into  a  tapd  loop  an  sho^m  in  71g.  2«3. 
The  mloroT^one  voltag©  was  flltorod  so  that  Its  speotnin 
was  that  passed  by  a  1/3  octavo  band  filter  v/lth  a  geo^ 
Taetrlo  aean  fr^ouoncy  of  4ooo  cpc.     Tbm  filtered  adcro- 
phono  volta^  Tias  fed  into  ohannel  one  of  the  oorrslator 
tine  delay  and  the  unfiltered  loudepeaker  voltaQs  was 
plaoed  acrons  channel  tvio  of  tho  correlator^     An  ViC  in- 
tegrator tiiae  of  1  second  and  a  tlae  delay  seanning  rat* 
of  1  oilHsocanc'  per  minute  v»r©  used*     The  resulting 
linear  plot  of  tho  cross  eorrelation  function  togothor 
with  the  -ilot  for  the  coraiKipable  pulce  tost  io  shoim  in 
Fig.  3«3.     In  "Wain  flr^ure,  oa<^  (sroup  of  pealEs  in  the 
corrolc^tion  functitm  ropresonts  the  arrival  of  the  sound 
from  an  ima^o  souroe  In  much  the  same  v/ay  as  tho  peaks 
in  the  accc8n:panylng  puis*  unalysis*     -An  accompanying 
^•8  l^lot  of  tho  correlation  v;as  made.     This  plot  is  not 
pflpee^nted  beimuse  of  the  extreme  physical  length  of  the 
plotting;  tape.     However,  Fig.  3»$  contains  a  pl0t  of  the 
peak  values  of  tho  eorrelation  curve  as  a  function  of    C 
f«r  Fig,  3»2«     Plotted  on  the  same  coordinates  is  a  graph 
of  tho  peak  pulse  paf^ssure  level  as  a  function  of  dlotanoe 
ti^VDled.     Oftie  points  for  "Uiis  seoond  curve  are  talcen  froa 
Fig.  4  (reference  8)» 
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irfTERSE   SQUARE   CITRVE 

EXPERIMENTAL   POINTS: 
A  -■     PULSE   METHOD 
O  -'-  CORRELATION  METHOD 


DISTANCE  TRAVELED   BY  SOUND   IN  FEET 
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Fls«  3*5     Plot  of  th©  tranfilent  reBpcnoo  of 

Room  10-390A  for  pulce  and  correlatlcfi 
analysis  * 
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3*4    DISCir-ION  OP  experhomtal  lasoLts 

Tho  data  on  pulae  statistics  aaaljals     sno^/s  the 
spaclnG  of  pialaee  arrlvlns  at  the  microphone  poeltlon 
vrtion  "both  aourco  andl  receiver  nr&  In  tho  comer  of  th» 
room*     The  correlation  analysis  of  tho  tranniont  respottM 
of  tho  room  eshould  show  eoraiparable  peaks  in  tho  correlation 
curTc  for  tliao  delays  eorre8i>onding  to  theno  i:n2lnG  nmelngs, 
!Ehe  ccrrGlation  curvo,   i-'ig.  3,3#  ahov/s  discroto  peolss  at 
delay  tlraos  of  12.7|  13»6,  21«2,   34.4,  37«2,  and  47.5  allll- 
seoondo*     rteforonco  to  Pig*  14  of  roferonee  8  efcowa  that 
these  delay  tiaies  corro  spend  roucshly  witri  tao  trancit 
tiaes  for  the  various  roflectod  pulses, 

Tb6  relative  amplltudQs  of  the  peaks  of  tho  corrola- 
ticai  Give  a  quantitative  aeaeure  of  th©  coabinod  offocts 
of  ( 

1.  At  tonus,  tl  on  of  sisnal  caui^d  by  InvorBO  sqxxare 
djjainuation  of  radiated  sound  intensity, 

2,  Attenuation  of  signal  caused  and  deereaao  la 
cohcrenee  of  radiation  fran  the  various  imase 
souroes  caused  by 

a)  Tho  transfer  function  of  one  rofloctlns 
sixrfaea  for  the  caae  of  ncm«de0m»]*at# 
dlBcrote  firct  order  laas*  ttourees. 

b)  Tbm  oeabined  transfer  function  of  all  the 
r© fleeting  surfaoes  for  higher  order  non- 
Ai^tnerate  Iraage  sources* 
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o)  Oior«ABe  in  ooherenoo  of  the  radiation  froa 
groups  of  degenerate  im'^^e  eotiroee. 
Th\is,  th»  plot  of  the  correlation  function  va  time  delay 
is  a  aeaauro  of  the  short  tem  transient  vu^pODBm  of  the 
room  for  a  parti ctilar  source  and  receiver  i^ooitlon. 

tHeins  correl3.tion  analysis,  it  is  also  ppeoiblo  to 
derive  a  value  for  reverberation  time  based  on  the  transient 
respcnre  of  the  room*     The  reverberation  time  siraply  eor- 
reepondc  to  the  difference  in  delay  tioo  bet^ieen  the  peak 
in  tho  curve  correspcjidins  to  the  arrival  of  sotsid  from 
eouroe  0»  0,   0  and  a  peak  in  tho  correlation  cui^vo  vrtoloh 
la  one-millianth  of ,  or  6o  db  dovm  frora  the  peak  of  the 
direct  ecomd,     l^is  time  may  be  found  by  oxtrapolr-ting  the 
rate  of  decay  indicated  by  the  short  term  transient  responne. 
For  the  example  discussed  in  the  prevlotie  section  \.1ier© 
both  oouroe  and  reeelver  were  in  the  cornor,  tlie  correla- 
tion curve  decayed  20  db  in  the  interval    ^  =:  10  ailll- 
seoondc  to^=  6?  mllllaeccjnds,     The  reverberation  time 
based  on  this  rate  of  decay  ie  3.1A  secorjds,     ^Hiis  value, 
then,  is  a  measure  of  the  revorberati on  tioe  for  a  1/3 
Mtave  bond  signal  having  a  geometric  mean  froquonoy 
«r  4000  cpa. 

In  line  i^ith  tho  experimental  pulce  analysis  of 
refOronce  8,  culditional  transient  reeponne  data  corres- 
pondlnG  to  the  data  of  Flss,  7  through  l6  has  boon  col- 
lected Cii^.  reduced  by  means  of  tlie  analog  correlator. 
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Tills  additional  datft  with  explanntory  notoa  lo  preflented 
in  Figs.  3.6,  3.7»  3.8,  and  5.9»     For  the  convenlenoa  of 
tho  reader,  a  copy  of  rof^wmc©  8  Is  Included  In  th» 
Appendix* 
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IV.    MXkstnaxans  or  dippusion 

4.1  INTRODUCnOH 

In  the  literature  on  room  acoustics,  much  attc  itlon 
hae  boon  directed  toward  defining  properties  of  ■^o\mCL 
fields  nhl<di  oeom  to  hr^ve  a  perceptible  and  qtiaatliatlro 
effect  upon  the  soodnesa  of  the  acoustical  envirosiznont . 
One  such  property  ie  •diffusion,*  without  dolvine  into 
the  tochnical,  psyoho-aooustic  aspects  of  tht  probKaaa, 
it  vould  \m   advisable  to  discuss  Bcme  of  the  more  intui- 
tively porceptibia  aspects  of  acoustics  which  are 
attendant  to  the  property  of  diffUBion. 

A  sound  field  is  defined  ao  being  ccsapletoly  diffuse 
If  it  has  uniform  energy  density  within  the  roGion  con- 
elderod,  and  if  the  direction  of  projxjsation  of  the  wave 
f rants  at  any  arbitrary  selected  points  are  wholly  randcaa 
in  diotributian*  In  discucsins  some  of  the  trends  in  the 
study  of  sound  waves  in  rocttt   SfexHie  and  Bolt  say,  *It 
now  appears  tliat  revert>eration  time  alone  is  not  always  a 
KUfficiont  aoasure  of  auditox^ium  excellonoo.  It  in 
desirable  to  h-rm   a  moan  squaa:^  presstu^  as  nearly  unifom 
as  possible  over  the  eeatllig  area*  It  is  also  important 
to  have  at  loaet  a  oez*tain  percentage  of  the  soxmd  reach 
the  listener  directly  from  the  speaig^r,  and  less  than  a 
certain  percentaoi  rea<^  the  listeners  indirectly,  after 
i:«flection  fraa  any  single  surface j  •..  "  In  the  lisht 
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ot  the  definition  of  diffttMntsSf  tht  tboiw  onotatioii 
px^soribus  that  tor  &r^  BwmA  field,  there  is  an  optlnnai 
degTM  of  diffttBiaii  liiioh  prodiioes  a  good  aoouetiool  en- 

YiMnmrnt.     Ons  may  reason  tliat  this  optimuB  decree  of 
diffusion  vd.ll  doT>«nd  upon  the  une  to  which  tiio  ooixnd  field 
is  to  bo  fxAopUA  i,e«  lecture  hall,  broadcast  otudlOi 
<laiittt\»  nmnio  hall|  and  more  particularly,  upon  tho 
••oastioal  tndltlon  of  tlie  population  in  the  sound  field. 
For  inntance,  the  population  of  Britain  may  be  said  to  have 
ill^tly  diffisrvnt  preJiwaee  ^haa  the  riopiidatioR  of  tlio 
Dhited   :tatea«     19)0  se  prefersaoes  would  probably  bo  based 
upon  tJie  way  the  British  have  been  uaed  to  hearing  their 
Biusic  or  locturee  and  also  upon  the  pec\aiaritioe  of  the 
talfttAfle*     the  degree  of  diffusion  in  any  Bound  field  will 
in  gcnora.1  bo  a  function  of  the  location  of  the  soui^ce,   tho 
MMipe  of  the  boundaries  enclosing  the  sound  fields  and  the 
shape,  iilM«99illett  eoeffident  and  distribution  at  the  siir- 
faees  in  the  aoond  field, 

h&  a  sort  of  intuitive  example  of  tho  effects  of 
diffusion  lot  us  vlsualir.e  a  full  symphony  orchestra  play- 
ing in  an  open  field.     Generally,  the  conductor  v;ould  be 
able  to  hear  clearly  only  those  instrtaaOTits  vAiiOh  v/ore 
near  to  him  or  %ihose  radiation  was  directive  and  boaaed 
in  hio  direction.     The  musicians  would  be  ablo  to  hear 
proportionally  less  of  the  total  effect  being  produced  by 
tho  ordiootra  as  a  group. 
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Under  thee*  clrcuoctanoos.   It  Ic  probable  tliat  a 
very  poor  porfonnanco  would  bo  the  rosxilt  of  this  lack 
of  dlffualon  or  nixing.     On  the  othor  hand,  If  oay  a 
MMlI  otriae  onsoiablo  tram  perfoxmlne  in  a  vory  diffuse 
biroftdoast  etiidio^  eaoh  aeniber  of  the  ensemble  would  be 
able  to  hoar  the  coabinod  effect  of  the  croup,  but  ha 
would  hear  rola timely  little  of  tha  soimd  that  hio  par- 
ticular inGtrument  >faa  producing  •     Suifii  a  conditicjn  would 
again  result  in  an  inferior  performance,      .'ith  the  proper 
degroo  of  diffuaion^  ho^foTer,  it  la  thoorizod  that  con- 
dltlone  v/ould  be  optioum  for  the  production  of  the  deelred 
aoouGtical  effect. 

In  a  technical  report  Number  B  058,   serial  No,  1953/ 
29*  dated  October  1953»  the  Brltl^  Broadcasting  Ccapangr 
discuG-ed  tho  property  of  diffusion  of  eound  fie Ida,     It 
vaa  said  that  isperfect  diffusicsi  aay  a;)pear  in  at  loaot 
tvro  moasu^rabXe  forms,    (1)  as  a  ohaiige  In  pros  euro  vrith  a 
chance  in  position,  the  frequency  of  the  ooimd  romalning 
ctMMJtant;    (2)  as  a  taiango  in  pressure  vrlth  change o  in  fro- 
qiiency,  the  positioais  of  source  and  rocolver  ronainins 
unchansod.     It  was  further  stated  that  imperfect  diffusion 
indicitod  in  the  steady  otn,te  aeasiirement  by  vrirtation  of 
pressure  vdth  frequency  or  pool tl on  is  shown  oy  aevlatlon 
in  tho  pressure,  time  rolationo  froa  an  cxi^onentlal  decay 
curve,     !Ehese  statements  are  perfectly  reasonable;  hovfever, 
the  functional  relationship  between  the  desree  of  diffusion 
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Urn  traaai«BLt  rBsponae  of  a  soOBd  field  is  not  ap> 
par»nt« 

Correlation  teoliniquM  B^&m.  to  off«r  a  method  of 
obtaining  a  aeooitzv  of  th«  dlfftiBlon  of  a  sound  field 
\f4hieh  depends  diroetly  upcn  the  definition  of  dlf f union • 
Wimnirrr,  this  method  esploya  almoet  the  cane  eqiaipaent 
•ad  armaiBBSionts  aa  are  used  for  the  tronoient  reepixise 
aeaaurementG  * 

4«2     COmn.ATKM  HSI^OD  07  DIFFUSION  mASimifmT 

ijot  un  Tisualize  the  following  experlaant, 
Me  ahall  place  a  microxiione  in  the  sowift  field  created 
by  a  loudGpealwr  driven  by  a  wide  bend  noise  voltafii 
generator* 

Now  let  tis  auto-eon^late  the  output  voXtag©  of  the 
microphone.     If  the  noise  voltage  has  an  infinitely  wi<^ 
band  flat  apoctrvai»  tiion  it  haa  been  ehovm^  that  tlio  a\xto- 
eorrelatian  ftfiieti<»i  of  the  !nicr<isliane  volta^a  ic  exactly 
the  combined  impulse  resp«ise  of  t^e  loudspealssr  end  the 
ndorophone.     ?*or  noise  voltajses  of  finite  band*.ddth,  this 
auto-correlation  fimctlon  is  the  transient  responr,©  of  the 
ejrstes  to  a  pulse  having  the  correfipondinc  spectral  coo- 
p<Wiltlcm.     In  ^neral,  thle  auto-oorrelaticm  function  haa 
*  WKKSmm  peftle  for  zero  time  delay  and  a  first  asdo  eroe* 
sins  ^^  a  tisje  delay  dependent  upon  the  bandwidth  cf  the 
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A01b«  voXtae^  (»•  -  CH  fluod  the  quantity  oo  -t*  a>i*     A 
••rl#s  of  Ruch  uuto-ooritilatiana  la  shcywn  In  ?ig«  A.l« 
As  Ax-olalnod  In  Chapter  1,   tha  bandtfldth  and  tho 
first  axis  cz*08slxig  of  tho  auto-oorrelation  ximotian  vary 
lanr»r«ely«     That  le^  tho  wider  the  band'^dth,   tho  ohortor 
Vm  dlffbrene*  In  d«lay  tinea  bot\n9on  tho  -neak  of  the  oor- 
rftlati on  fuactlon  ajoA  Urn  first  axis  oroQalng*     It  It 
pfMslble  to  meaoitre  the  tlae  delay  corrospoodlng  to  this 
first  axle  crogalng  with  a  hlgli  dmigpme  of  acovuraey.     Thm 
tino  aolay  roglster  assoolated  with  tho  analog  corrola- 
tlon  oonputor  po  units  roadlnga  of  time  delay  to  the  tenth 
•f  a  mllllnocond ,     It  Is  posolble  to  Interpolate  b«t«ift«R 
Urnue  rogl^ter  readings  by  the  follo^^lag  aotliod, 

a)  Into  eaoh  channel  of  the  tlste  delays  feed  a 
•tiiadard  10  Scoxm  signal  voltage. 

b)  mmBVgt^  the  jjkuum  dlff^renoe  of  Hie  outputs  cat 
the  tlae  dolay  for  thle  10  kops  signal, 

c)  Slne^  the  period  of  the  10  keps  signal  is  pre- 
cisely l/lO  of  a  mllllseeond,  it  is  pooclblo  to 
hand  adJTist  the  tliae  delay  siefihanlsm  to  a  valtid 
ewr^spoadlng  to  any  fractlcai  of  a  period  In 
terns  of  diigroos  of  phase  difference  beti^eon  th» 
tla©  delay  outputs, 

d)  ^^%h  the  tliae  delay  meohanlsffi  thus  adjusted,   csub 
Oaa  amr  remove  the  10  kcpo  signal,  apply  tlie  nois^ 
voltage,  and  meaGure  the  voltage  output  correapood- 
Ing  to  the  value  of  the  auto-oorrelatlon  curve  for 
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that  tiae  aolay. 

RoW|  knowing  th»  tUm   d^l&y  of  the  axis  orossing, 
one  can  plaoe  two  such  identical  mlcropiiocioe  In  thA 
sound  field  of  a  souroe  radiating  epherloal  soond  waTes. 
If  the  mlcrorihaiee  are  lilaeed  a  dlstsnee  apart  equal  to 
the  tlae  delay  to  the  first  axis  ez*o88lns  of  the  auto- 
eorrelatlcn  curve,  and  then  orientated  In  such  a  way  that 
the  sound  fronts  frca  the  source  strike  both  milDDs  at  the 
eaae  Instant  (Pig*  4*2) |  then  the  cross -eorrelatl en 
fimctlcn  at  zqvo   tlae  delay  for  the  outp^  Toltefls  of 
the  two  alcrophoces  will  correspond  to  the  value  of  the 
auto-correlatlon  cimre  for  «»  0.  However,  If  the  couroe 
mlcrophoaae  array  Is  turaed  90  degrees  In  a  horlsontal 
plane  so  that  the  sound  wares  strike  one  microphone.  Is 
subjected  to  a  physical  delay  corresponding  to  the  axis 
crossing  tlae»  and  then  strlkos  the  secoood  microphone 
(Fig*  ^.3)»  then  the  corresponding  cxvjss-correlatlott 
function  will  be  zero.  Thus  It  Is  poeclble  to  obtain  a 
plot  of  cross-correlation  vs  l^e  orientation  the  nlcro- 
fticnsc  array  in  degrees  from  0  to  90° «  If  we  specify 
that  sero  degrees  eorrespends  to  the  orientation  vdien 
both  sites  are  equidistant  from  the  source,  then  we  oan 
•flflrert  our  preriouely  meaaxired  aecurate  plot  of  the 
auto-eorrelation  cutto  tb  Umb  delay  to  a  plot  of  cross- 
correlation  function  vs  orientation  In  degrees  by  means 
of  the  f  oznula  =  d  sin  G  where  d  Is  the  milDe 
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fie*  4*2      ItMrtograph  of  «lcrophoi»  array  oriented 

at  0  a06f««»  In  smma  field  of  loii^ispoater. 
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Us*  4*2      aiiytograpli  ot  adcrophono  fovo^  oriented 

at  0  assro««  In  BcnmA  field  of  loudspoafaar. 
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Big,  4,3       aiotograTii  of  nderophooe  array  oriented 

at  90  degrees  in  sound  field  of  loudspeaker, 
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separation  In  feet,   c  is  the  speed  of  sound  In  feet  per 
millisecond,   and       la  time  delay  In  milliseconds.     These 
calctLLatlons  have  been  plotted  for  a  nolne  spectrum  flat  from 
100  cps  to  10  kcps,  having  a  roll  off  of  6  db  per  octrv^^ 
beyond  these  frequencies   (Fig.   4.4). 

The  microphone  array  together  with  the  multiplying 
and  integrating  circuits  of  the  analog  correlator  comprise 
an  instrument  Vihlch  gives  a  measure  of  the  properties  In- 
volved In  the  definition  of  diffusion. 

Suppose  the  microphone  array  v/as  placed  in  a  perfectly 
diffuse  field.     Then,  no  matter  which  way  the  array  v/as 
t\imed,   ttie  value   of  the  cross-correlaticai  function  woiild 
be  the  same,   since  by  definition,   a  diffuse  field  is  one 
in  which  th©  wave  fronts  appear  to  be  coming  uniformly 
from  all  directions.     It  is  clear  that  for  variations  in 
diffusion,   the  related  cross-correlation  CTirv3  vrould  have 
some  chr.racterlstic  shape. 

4.3  EXPERIMEKTAL  RESULTS 

EXPERIKSKT  1 

T^'7o  illtec  21-3R-150A  microphones  vrere  placed  two 
inches  apart  in  anecholc  space  in  the  sound  field  pro- 
duced by  an  Altec  728-B  closed -box-baffled  loudspeaker. 
The  cross  correlation  of  the  clipped  output  voltage  of 
the  microphones  was  plotted  as  a  function  of  orientation, 
Fig,  4,5,  from  0  degrees.  Fig.  4.2,  to  90  degrees, 
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yig.   ^.3«     Thei   curve  In  Pic  4.5  Is  a  meanur©  of  the  dif- 
fusion of  anecholc  spa^ce.     As   should  be  expectod,    the 
source  shov;s  that  the  sound  field  produced  by  the  loud- 
speatosr  In  anecholc  space  Is  not  diffuse. 

The  value  of  microphone   separation  uced  In  this  first 
experiment  v/as  obtained  by  trial  and  error.     It  Is  impos- 
sible to  get  the  desired  results  by  spacing  the  micro- 
phones the  theoretically  correct  distance  apart  because 
of  the  finite  size  of  the  microphone  diaphragms. 

EXPERIMENT  2 

The  following  experimental  work  was  performed  in  the 
model  stiadlo,  Rooa  20-F-009A,  located  in  the  Acoustics 
Laboratory  at  M.I.T.     This  studio  was  designed  as  a 
diffuse  room.     A  physical  description  of  this  room  is 
contained  in  reference  12. 

With  the  loudspealoar  facing  the  v/all  shovjn  in  Figs. 
4,2  and  4.3#  diffusion  data  was  taken  with  the  microphone 
array  placed  at  various  distances  for  a  line  perpendicular 
to  the  radiating  face  of  the  loudspealssr  enclosure.     The 
noise  voltage  and  other  equipment  v/as  the  same  as  in 
Experiment  1.     As  in  Experiment  1,   the  microphone  output 
voltagos  were  clipped  20  db  below  their  rms  value. 

The  results   of  the  diffusion  measurements  together 
^vlth  the  plot  of  the  transient  response  of  the  studio  are 
sho^-m  in  Figs.   4.6  and  4.7  respectively. 
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EXPERIMENT  3 

The  follov/lng  experimental  data  was  tal^n  In  Room 
10-390-A,   the  hard  plaster  roam.     In  these  experiments 
the  data  v/as  tape  recorded  and  then  "brought  into  the 
laboratory  and  reduced. 

The  same  equlpjaent  used  in  Exiperlmento  1  and  2  was 
used  in  the  follo^dng  tests. 

Diffusion  data  and  tranoiont  response  data  \tb.s,  taken 
for  various   ccnditione   of  the  room.     For  these  terts,   tlie 
loudspeaker  v/as  placed  facing  a  comer  of  the  rocra  and  the 
diffusion  data  ^^^a3  taken  for  various  positions  in  the  room. 
The  positions  were  located  at  distances  of  20,   15,   10,   8, 
and  5  feet  from  the  loudspeaker  on  a  line  bet^^reen  the 
loudspeaker  location  ajtid  the  diagonally  opposite  comer 
of  the  room.     The  transient  response  data  vxas  taken  at  the 
15  foot  position. 

The  above  described  data  vms  taken,   reduced,   and 
plotted  for  the  following  rocaa  oonditionG:  ' 

a)  All  walls  of  the  rocm  were  bare.     Figs.  4,8  and 
4.9. 

b)  One  long  vraill  of  the  room  was  covered  id.th  a 
highly  absorbent  Plberglas  filled  quilt.     Pigs. 
4.10  and  4.11. 

c)  One  long  vrall  and  one  sliort  v/all  v/ere  covored  \ilth 
quilting.     Figs.   4.12  and  4.13. 

d)  T%-70  short  walls  and  one  long  wall  v;ere  covered  vrlth 
quilting.     Figs.  4,14  and  4.15. 
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4.4  DISCUSSION  OF  SXPEriDIIilNTAL  RESULTS 

In  reviewing  the  experimental  results  of  3eotlon 
4*3  it  appears  that  the  measuring  system  and  the  tech- 
niques enployed  yield  results  similar  to  those  obtainable 
v:lth  a  hishly  directiona.l  microphone.     It  appears  that 
the  practice  of  setting  the  microphones  a  cer*tain  dis- 
tance apart  corresponds  to  maVring  a  microphone  which  is 
extremely  serjsitive  in  measuring  the  direction  of  prop- 
ogation  for  a  sound  having  a  particular  siDeciilar 
campoGltion*     Generally,   setting  the  microi^ones  further 
ajjart  maJaes  the  array  more  sensitive  to  spectra  having 
lower  Geometric  mean  frequencies.     Insxifficient  data  is 
available  to  permit  an  evaluation  of  this  method  of 
measurins  diffusion.     However,   the  author  believes  that 
the  data  presented  herewith  is  sufficient  to  indicate  a 
possible  experimental  method  of  obtaining  a  laeasxiremont 
of  the  degree  of  diffusion  of  a  sound  field.     Moreover, 
it  is  believed  that  this  method  i/arrantc  fTiirbher  study, 

4.5  StJQGESTICffilS  FOR  ADDITIONAL  V/ORK 

On  the  basis  of  the  discussicaa  and  measurements  of 
diffusion  made  in  cannection  v/lth  this  project,  it  is 
felt  that  there  is  some  potential  merit  in  the  techniques 
employed.  As  a  suggestion  for  addlticnal  v/ork,  it  Is  pro- 
posed that  some  consideration  be  given  to  further  diffusion 
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measurement n  made  x/lth  varloun  "bnndxTl.dthr,   of  noise  voltage, 
Such  a  projGct  v;ould  mafce  errtenslve  experiment g  in  some 
roan  such  as  the  model   studio.     One  possible  series  of 
experiments  would  Involve  mea-suremente  for  various  micro- 
phone  separations   (selected  In  a  manner  similar  to  that 
of  Experiment  1,    Section  A. 3)  aJid  various  conditlcns  of 
diffusion  (adjustiaents  in  diffusion  to  be  made  by 
rearran^^lng  variable  diff users  in  the  model  studio) , 
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V.      CONCLUSION 

5.1     SIB^MAiii  UF  CORRELATION  MEASUREMENTS 

CoiT^Ql'-iticm.  measuremoiits  as  studied  in  this  project 
have  been  of  two  types   -   (1)  msasiirements  of  tho  transient 
respoiiGe  of  rooms  and   (2)  measurements  based  on  the 
defiioitlon  of  a  diffuse  sound  field.     It  Is  an  interest- 
ins  coincidence  that  both  types  of  measurements  can  be 
made  vrltli  practically  the  same  experimental  setup. 
Transient  resi^onse  data  may  be  taken  by  cx»os8  correlating 
betvflsen  tho  microphone  output  voltages  of  tv/o  microphones, 
one  of  vjhlch  is  located  near  the  source  and  the   other  in 
the  reverberant  field  of  the  room.     Diffusion  measurements 
are  mad©  in  iiie  same  way,  by  cross  correlating  bet'^/ecn  the 
microphone  output  voltages  of  tv/o  microphones  which  are 
spaced  a  particiilar  distance  apart  and  oriented  in  a 
certain  v/ay  with  respect  to  the  source* 

The  above  illustrated  coincidence  serves  to  point  out 
the  versatility  of  correlation  tedrinlques.     Unfortunately, 
this  coincidence  also  serves  to  illustrate  one  of  a  number 
of  pitfalls  that  one  may  encounter  In  the  use  of  such 
techniques.     Improper  microphcaae  placement  in  either  type 
of  measurement  would  lead  to  spurious  results. 

It  was  stated  before  that  diffusion  meacurements  by 
correlation  technlqxjBs  seemed  to  yield  results  similar  to 
those  obtainable  by  use  of  a  highly  directional  microphone. 
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It  should  be  of  Interest  to  point  out  th'^t  spaclnG 
mlerophones  a  dlotancse  apart  which  is  calculated  from 
auto-correlatlon  data  corix)3pondc  to  calibrating  the 
micro:Thone  array  for  maximum  directionality  for  the 
particular  frequency  spectrum  of  noioe  voltago  being 
used , 

This  directive  microphono  array  constructed  by  cor- 
relation techniqijes  appears  to  have  two  major  lobes  in 
its  directivity  pattern.  'These  major  lobes  are  180 
degrees  apart.  The  calibration  process  concists  of 
spacing  the  microphones  so  as  to  eliminate  any  minor 
lobes  in  the  directivity  chej^acterlGtics  of  the  array. 

It  should  be  pointed  out  that  the  directivity  1*1! oh 
we  have  been  talfclng  about  refers  only  to  the  ability  of 
the  microphone  array  and  the  associated  correlation  equip- 
ment to  discriminate  betv;Ben  the  presence  of  a  plane  wave 
front  vrhlcii  reaches  both  microphones  at  the  same  Instant 
and  a  plane  v/ave  front  which  reaches  each  microphone  at  a 
different  instant. 

In  a  previously  mentioned  report  by  the  British 
Broadcasting  Company,  It  vras  concluded  on  the  basis  of 
their  experimental  studies  that  the  methods  of  short  pulse 
analysis  are  not  satisfactory  for  the  investigation  of 
diffusion  in  full  scale  rooms.  The  pulse  methods  refer- 
red to  in  tills  conclusion  involved  -  (a)  studies  of  the 
irregularity  of  the  envelope  of  the  decay  for  a  short 
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pulse  of  Boimd  and  (b)  variation  In  those  short  pulse 

decay  irregularities  \^th  position  in  the  room. 

l"? 

Another  group  of  inventlgfitors  -^  uain^  a  direc- 
tional micro^:ihone  technique  attempted  to  obt-i^-ln  a  cteady 
state  raoacuix)  of  the  diffusenoss  of  sound  fiexus.  The 
conclusion  to  their  researcli  cited  the  nocesslty  for  a 
microphone  that  could  be  made  hichly  directional  over 
the  whole  audio  f2?equency  range. 

It  appears  that  the  correlation  method  of  laeasuring 
transient  response  is  suitable  for  use  in  full  scale 
rooms,  Unlilse  short  pulse  methods,  correlation  analysis 
of  transient  response  may  bo  made  to  yield  short  tein 
decay  rates  for  various  bandv/idths  of  frequencies.  On 
the  basic  of  this  latter  fact.  It  may  be  concluded  that 
correlation  analysis  presents  a  connecting  linl^  botv.'een 
room  response  data  gathered  by  steady  state  methods  and 
tfbat  gathered  by  ptilse  methods. 

In  addition,  it  can  be  seen  that  correlation  tech- 
niques, v;h$n  u.^ed  in  conjunction  vdth  a  microphone  array 
calibrated  for  the  purpose,  can  be  made  to  yield  data 
similar  to  that  obtained  by  directional  micro]±.onB 
methods.  Further,  the  microphone  array  can  be  calibrated 
for  amy  bandi-rldth  by  use  of  data  talsBn  from  Its  auto- 
correlation curve. 

Pron  tiiese  facts  it  may  be  further  concluded  that  the 
correlation  method  of  analysis  of  bounded  sound  fields  is 
potentially  more  Informative  than  either  steady  state 
methods  or  pulse  methods  taken  alone. 
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Many  sounds  of  si)ecch  and  music  more  nearly  rcsemlile  pulsed 
wave  trains  than  aliruptly  terminated  continuous  sounds  as  used 
in  reverberation  measurement.  It  is  therefore  not  surjjrising  to 
find  that  two  rooms  can  differ  markedly  in  acoustical  riualit\' 
even  if  they  appear  identical  under  reverberation  analysis  which 
ignores  details  of  short  transients. 

This  paper  introduces  a  pulse  statistics  point  of  view  which 
takes  immediate  account  of  the  pulse-like  nature  of  common 
sounds.  Fundamentally,  the  method  consists  in  examining  the 
response  of  the  room  to  a  short  pulse.  The  walls  are  replaced 
by  an  array  of  image  sources  (simple  images  if  the  walls  are  hard, 
or  appropriately  modified  if  there  is  absorption).  These  image 
arrays  are  then  considered  statistically. 

From  this  approach  one  can  derive  such  classical  quantities 
as  reverberation  time  and  mean  free  path.  One  can  also  analyze 


the  detailed  nature  of  discrete  reflections  including  interference 
effects,  and  thus  obtain  an  average  correlation  between  room 
geometry  and  the  character  of  its  jjulse  response. 

Itlealized  cxjieriments  in  a  hard-walled  rectangular  room  are 
employed  to  illustrate  the  essential  features  of  this  approach. 
.•\  point  source  emits  an  exponential  damped  3600-c.p.s.  wave 
train  of  about  2  msec,  duration.  The  received  signals  are  recorded 
logarithmically  on  an  oscillograph  and  the  system  is  calibrated 
for  quantitative  results.  Several  dozen  discrete  reflections  can  be 
measured  and  correlated  with  calculation.  The  pulses  merge  into 
a  more  or  less  continuous  background  after  a  time  that  is  calcu- 
lated and  confirmed  experimentally.  Detailed  differences  arise 
according  to  the  positions  of  the  source  and  microphone  in  the 
room. 


INTRODUCTION 

TT  is  well  known  that  reverberation  time  is  not  a 
^  completely  adequate  index  of  the  acoustical  quality 
of  a  room.  Implicit  in  the  concept  of  reverberation  time 
is  the  assumption  that  the  room  reaches  a  steady-state 
diffuse  condition  before  the  source  of  sound  is  abruptly 
terminated.  In  practice,  however,  most  sounds  of 
speech  and  music  can  be  generally  classified  as  pulsed 
wave  trains  whose  amplitudes  and  frequency  compo- 
nents fluctuate  sufficiently  within  time  inter\-als 
shorter  than  the  time  constant  of  the  room,  so  that  the 
room  seldom  reaches  steady  state.  Thus  it  would  seem 
that  the  response  of  the  room  to  transient  sounds  of  this 
general  type  is  an  especially  important  physical  prob- 
lem of  room  acoustics.  The  results  obtained  recently 
by  Mason  and  Moir^  and  others-  who  have  used  short 
tone  bursts  to  investigate  acoustics  of  auditoriums  lend 
support  to  this  point  of  view. 

The  task  of  describing  mathematically  the  response 
of  a  room  to  an  arbitrary  transient,  and  of  studying 
the  roles  of  room  geometry  and  distribution  of  absorbing 
materials  in  this  response,  is  extremely  complicated. 
The  problem  can  be  approached  from  a  normal  mode 
point  of  view,'  or  one  can  attempt  to  "follow"  the  sound 
waves  around  in  the  room  as  they  are  reflected  back 
and  forth  from  the  walls.  The  latter  approach  has  been 
recently  investigated  by  Mintzer"*  using  Laplace  trans- 
form methods. 

From  the  transient  point  of  view,  it  is  desirable  to 


*  This  work  was  supported  in  part  by  the  ONR,  Department 
of  Navy,  under  Contract  NObs  25391,  Task  7. 

'  C.  A.  Mason  and  J.  Moir,  "Acoustics  of  cinema  auditoria," 
J.  Elec.  Eng.  88,  Part  III,  No.  3  (September,  1941). 

^  British  Broadcasting  Corporation,  Engineering  Division, 
Research  Department  Reports  B.027  and  B.035. 

^  P.  M.  Morse  and  R.  H.  Bolt,  "Sound  waves  in  rooms,"  Rev. 
Mod.  Phys.  16,  117  (1944). 

*  D.  Mintzer,  "Transient  sounds  in  rooms,"  J.  Acous.  Soc.  Am. 
22,  341  (1950). 


use  the  second  approach.  In  essence,  this  method  con- 
sists of  replacing  the  effect  of  the  boundaries  of  the 
room  by  an  infinite  array  of  image  sources,  each  image 
corresponding  to  one  of  the  multiple  reflections  of  the 
original  wave  emitted  by  the  source.  Finding  these 
images  analytically  is  no  simple  matter,  and  only  in 
ver>^  special  cases  will  the  images  be  "mirror  images"  of 
the  source.'  The  image  concept  has  been  used  exten- 
sivel}-  in  earlier  geometric  studies^' ^  where  the  source 
is  considered  to  be  incoherent.  Due  mainly  to  mathe- 
matical difficulties,  little  use  of  images  has  been  made 
as  yet  in  wave-acoustical  investigations  of  rooms. 

However,  it  has  become  increasingly  apparent  in 
recent  years  that  the  first  20  db  of  decay  of  a  sound  in  a 
room  is  of  primary'  importance  in  differentiating  be- 
tween two  rooms  which  have  approximately  equal  over- 
all reverberation  times.  Further,  as  the  work  of  Mason 
and  Moir  indicates,^  the  time  and  amplitude  distribu- 
tions of  reflected  tone  bursts  can  be  correlated  with  the 
acoustical  quality  of  a  room.  These  facts  indicate  that 
the  images  relatively  close  to  the  source,  i.e.,  the  first 
few  reflections,  are  primarily  responsible  for  certain 
important  features  of  the  acoustical  character  of  rooms, 
as  Brillouin  has  observed.^  It  should  therefore  be  worth 
while  to  study  this  "short  term"  transient  response  by  a 
method  of  images  in  which  all  wave  properties  of  the 
image  sources  can  be  considered  (i.e.,  where  the  assump- 
tion of  an  incoherent  source  is  not  made).  Further,  if 
an  image  array  satisfying  the  boundary  conditions  can 
be  found,  one  should  be  able  to  treat  this  array  sta- 
tistically and  thus  obtain  the  long  term  average  tran- 
sient response  as  well. 

This  paper  is  confined  for  the  most  pajt  to  a  discus- 
sion of  an  idealized  case,  a  hard-walled  rectangular 


'  C.  F.  Eyring,  "Reverberation  time  in  'dead'  rooms,"  J.  Acous. 
Soc.  Am.  1,  217-241  (1930). 

^  J.  Brillouin,  "Sur  I'acoustique  des  salles,"  Rev.  d'Acoustique 
1  (September-November,  1932). 


328 


329 


PULSE     STATISTICS 


room  containing  a  simple  source  which  emits  a  short 
pulse.  Statistical  properties  of  the  image  array  are 
investigated  to  illustrate  the  image  method  and  to 
provide  a  basis  for  future  work  on  more  general  cases. 
The  Umited  applicability  of  results  based  on  a  simple 
mirror  image  picture  has  already  been  pointed  out.^ 
However,  as  the  discussions  of  both  Morse  and  Bolt^ 
and  Mintzer^  indicate,  a  specular  approximation  is 
allowable  in  cases  where  the  walls  are  not  too  soft,  and, 
in  any  case,  one  can  assume  that  the  image  of  a  simple 
source  can  be  described  analytically  by  an  expansion  in 
spherical  harmonics,  which  is  essentially  a  multipole 
expansion  of  the  image  source  taken  about  the  mirror 
image  point.  Thus,  if  the  walls  are  not  fairly  hard, 
neither  the  simple  mirror  image  nor  the  specular  ap- 
proximation can  be  used,  and  other  suitable  representa- 
tions of  the  images  must  be  found. 

PULSE  STATISTICS  THEORY 

1.  Image  Space  and  the  Time  Distribution  of 
Reflected  Pulses 

With  these  restrictions  in  mind,  let  us  now  set  up  a 
working  picture  for  pulse  analysis.  We  shall  consider  a 
simple  rectangular  room  with  perfectly  reflecting  walls. 
The  room  has  dimensions  Lx,  Ly,  L^,  one  comer  being 
at  the  origin  of  a  Cartesian  coordinate  system. 

A  sharp  pulse  of  sound  is  emitted  from  a  source  in 
this  room.  The  source  is  assumed  to  be  a  simple  one, 
with  spherically  symmetric  radiation.  This  point 
source  may  be  located  by  the  vector: 


r,  =  Xi+Fj+Zk, 


(1) 


is  located  at  a  position : 


rR=U\+V]+Wk. 


(2) 


The  vector  displacement  of  the  receiver  from  the 
source  is : 

Ro=r,-r/e.  (3) 

Associated  with  this  sound  source  is  an  infinite  array  of 
image  sources  each  occupying,  alone,  an  image  replica 
of  the  physical  room.  This  image  array  also  is  illustrated 
in  Fig.  1.  Each  image  cell  is  designated  by  three  numbers 
(/,  m,  n),  and  these  three  numbers  can  take  all  integral 
values  from  minus  infinity  to  plus  infinity.  The  cell 
(0,  0,  0)  is  the  actual  room. 

The  pulse  associated  with  the  image  (1,  0,  0)  reflects 
once  from  the  wall  x=Lx.  The  pulse  from  (2,  0,  0)  re- 
flects first  from  the  wall  x=0,  then  from  the  wall  x=Li. 
Obviously  the  absolute  value  of  the  cell  number,  /,  for 
cells  lying  along  the  x  axis  gives  directly  the  number  of 
wall  reflections  suffered  by  the  pulse  from  the  image  in 
question.  The  pulse  from  (1,  1,  0)  reflects  once  from  the 
wall  x=Lz  and  once  from  the  wall  y  =  Ly,  so  that  its 
number  of  reflections  is  |  / 1  +  |  w  ] .  In  fact,  by  this 
system  of  designation,  the  total  number  of  wall  reflec- 
tions suffered  by  the  pulse  from  the  image  (/,  m,  n)  is 
directly : 


N 


Imn ' 


U   +    >?^|+  1^1 


(4) 


The  vector  position  of  each  image  source  is : 

rimn  =  x/i+y,„j+s„k.  (5) 

The  vector  position  of  each  image  with  respect  to  the 
receiver  is : 

R:m7i=  ffmn— Tie.  (6) 


as  shown  in  Fig.  1.  A  point  receiver  of  sound  pressure      To  evaluate  these  last  two  equations,  we  note  first  that 


Fig.  1.  An  image  ar- 
ray, showing  image  cells 
and  source  and  receiver 
locations.  The  cell  desig- 
nation numbers  are  en- 
closed in  the  small 
boxes.  (Z,x  =  23.0  ft., 
Z,„=13.4  ft.,  L,=SM 
ft.) 
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each  image  cell  is  a  mirror  reflection  of  the  cells  ad- 
jacent to  it.  This  complicates  the  analysis  somewhat 
in  that  we  have  two  kinds  of  symmetry  with  respect  to 
the  basic  coordinate  system.  Fortunately  the  cell 
designation  numbers  indicate  directly  the  type  of 
symmetry:  even  numbers,  including  0  for  the  actual 
room,  designate  cells  in  which  the  source  position 
duplicates  that  in  the  original  room,  while  odd  num- 
bered cells  have  sources  at  a  reflected  position.  There- 
fore the  components  for  Eq.  (6)  are  given  by : 


symmetry-.  From  Eq.  (6)  we  then  get: 


xi=iLr-\-x  ^ 


xi=(l  +  l)L.-X 


-y\ 


ym={l+m)Ly 
Zn={\+n)Lz—Z 


I,  m,  n  even, 


/,  m,  n  odd. 


(7) 


In  order  to  generalize  some  illustrative  calculations  we 
select  the  longest  dimension  of  the  room  as  a  scale 
unit  and  introduce  tlie  following  definitions : 


Lx=L,     Ly=pL,  Lz  =  qL,     p,q<l; 

(£/-X)/Z,=  Mxe,  (f/+X)/L  =  Mx<„ 

(F-F)/i:=Mv«,  (F+F)/Z  =  M„o, 

(PF-Z)/I  =  Mze,  (IF+Z)/L=M.<,. 


(8) 


Thus  p  and  q  are  dimension  ratios  and  the  ju's  give  the 
source-to-receiver    displacements    for    both    kinds    of 


R/mn 

L 

[/-/i„]i 
or 

>  +  ■ 

or 

• 

[«V-MrJk      ■ 

even 

+  ■ 

or 

-,      ■ 

or 

[(nn 

hl)<7-M.o]k, 

odd 

(9) 


Some  calculations  from  this  equation  are  illustrated  in 
Figs.  2  and  3  which  will  be  discussed  later. 

Next  let  us  find  the  average  number  of  pulses  re- 
ceived up  to  a  specified  time,  /,  after  the  emission  of  the 
pulse  from  the  original  source.  These  pulses  come  from 
all  of  the  images  within  a  radius  |  RimnI  =d.  One  image 
is  contained  in  each  cell  of  volume  V=LiLyL^.  The 
number  of  pulses,  A'p,  is  directly  given  by  the  volume 
of  the  sphere  out  to  cl  divided  by  the  volume  of  one  cell : 


A^ 


47r|R,„„p     ^T(?^ 


W  3F 

and  the  number  of  pulses  received  per  second  is : 
dNr,    47rc^/2 
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Fig.  2.  Graph  showing  N p,  the  number  of  pulses  arriving  at  the  receiver  up  to  a  time  /  as  a  function  of  the  dimen- 
sionless  time  parameter  r.  Note  the  extreme  stepped  behavior  of  the  dashed  curve  for  a  cubical  room. 
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These  equations,  (10)  and  (11),  are  valid  for  values  of 
/,  m,  and  n,  large  enough  so  that  the  details  regarding 
positions  of  source  and  receiver  in  the  room  can  be 
neglected.  These  equations  also  are  illustrated  in  Figs. 
2  and  3.  A  calculation  of  an  actual  case  illustrates  that 
a  surprisingly  large  number  of  reflections  per  second  are 
predicted.  In  a  room  of  10,000  cu.  ft.,  Eq.  (11)  gives  a 
rate  of  180  pulses  per  second  at  1/100  of  a  second. 
There  are  several  conditions  met  in  practice  which 
greatly  reduce  this  number.  For  one  thing,  the  floor 
is  generally  quite  absorptive  when  an  audience  is 
present,  so  that  half  of  the  spherical  volume  just 
assumed  is  essentially  eliminated.  Also  there  are 
usually  a  large  number  of  degeneracies  or  coincidences 
which  reduce  the  effective  number  of  pulses,  as  we 
shall  see. 

Figures  2  and  3  show  yVp,  the  number  of  pulses  re- 
ceived up  to  a  time  /,  for  rooms  of  several  dimension 
ratios  and  various  positions  of  source  and  receiver. 
Np  is  plotted  for  convenience  as  a  function  of  the  dimen- 
less  time  parameter 


d 


ipq)iL 


(12) 


where  p,  q,  and  L  are  defined  in  Eq.  (8).  The  heavy 
solid  line  in  both  figures  gives  the  average  value  of  N p 
as  computed  from  Eq.  (10). 

In  Fig.  2  the  dashed  curve  shows  the  extreme  de- 
generacy of  a  cubical  room  when  both  source  and  re- 
ceiver are  in  a  corner.  The  pulses  arrive  in  large  groups 
because  of  the  symmetry  of  the  image  point  lattice  for  a 


cubical  room  and  the  location  of  source  and  receiver. 
The  dash-dot  curve  for  a  rectangular  room  with  /»  =  0.5 
and  9  =  0.33  shows  similar  but  smaller  groups  of  coin- 
cident pulses.  Both  these  curves  start  at  iVp=8  since 
the  direct  pulse  and  the  pulses  from  the  seven  source- 
corner  images  reach  the  receiver  simultaneously  at 
/=  r  =  0.  The  light  solid  curve  is  for  the  same  rectangular 
room,  source  and  receiver  now  being  at  point  such  that 
coincidences  are  more  or  less  "accidental."  We  note 
that  this  curve  follows  the  average  curve  very  closely 
for  T>2.2.  For  r<2.2,  the  curve  lies  consistently  below 
the  average  curve.  This  is  partly  because  the  separation 
of  source  and  receiver  are  such  that  no  pulse  arrives  at 
the  receiver  until  t=1.1.  In  detail,  this  initial  part  of 
the  curve  depends  in  a  somewhat  more  involved  way 
on  the  relative  positions  of  the  real  and  image  sources 
and  the  receiver. 

Figure  3  shows  in  more  detail  the  initial  rise  of  the 
curve  just  discussed.  The  dash -dot  and  dashed  curves 
show  similar  initial  rises  for  a  room  of  />  =  0.8  and  9= 0.6. 
The  same  general  behavior  is  shown  by  these  curves, 
although  both  reach  the  average  curve  sooner.  The 
duration  of  this  initial  departure  appears  to  decrease  as 
the  room  becomes  more  nearly  cubical. 

2.  Some  Special  Cases  of  Coincidences 

In  Eqs.  (10)  and  (11)  it  was  assumed  that  all 
images  are  distinct,  that  is,  that  no  two  pulses  arrive 
at  the  receiver  at  exactly  the  same  time.  However,  if 
source  and  receiver  are  both  placed  in  certain  locations, 
it  is  possible  for  the  pulses  to  arrive  in  groups,  as  is 
evident  from  Figs.  2  and  3.  It  is  of  interest  to  find  the 
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Fig.  3.  Graph  on 
an  expanded  scale 
showing  the  initial 
behavior  of  Np  as  a 
function  of  t. 
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number  of  distinct  pulse  groups  reaching  the  receiver 
for  some  of  these  degenerate  cases.  In  all  cases  it  is 
assumed  that  the  room  is  rectangular,  hard-walled, 
and  of  incommensurate  dimensions. 

Case  1 

Consider  the  source  to  be  in  a  corner  and  the  receiver 
in  a  position  incommensurate  with  respect  to  the  images 
but  otherwise  arbitrary.  If  the  source  is  close  enough 
to  the  corner,  the  images  will  clump  together  in  groups 
of  eight,  being  so  close  together  that,  efifectively,  each 
clump  of  eight  acts  as  a  single  distinct  image.  Since 
pulses  reaching  the  receiver  will  arrive  in  groups  of 
eight,  the  elective  number  N p  of  distinct  pulses  arriving 
per  second  is : 

1  Aiiicty    wicty 

(13) 


^y;= 


8     3V 


Case  2 


6V 


Here  we  place  both  source  and  receiver  in  the  center 
of  the  room.  All  the  images  are  distinct,  but  because  of 
the  symmetry  arising  from  the  position  of  the  source 
and  receiver,  a  pulse  from  an  image  along  a  negative 
coordinate  axis  (the  origin  being  taken  at  the  center  of 
the  room  for  convenience)  will  arrive  simultaneously 
with  a  pulse  from  an  image  along  the  positive  coordi- 
nate axis.  A  pulse  from  a  non-axial  image  in  the  first 
quadrant  of  a  coordinate  plane  will  arrive  simultane- 
ously with  pulses  from  corresponding  images  in  each 
of  the  other  quadrants,  and  a  pulse  from  an  "oblique" 
image  will  arrive  simultaneously  with  jjulses  from  cor- 
responding images  in  the  other  seven  octants.  Tlius  the 
effective  distinct  images  are  all  contained  in  one  octant 
of  image  space.  Counting  of  the  number  of  distinct 
pulses  arriving  can  then  be  accomplished  in  exactly 
the  same  manner  as  the  counting  of  the  number  of 
normal  modes  with  frequencies  less  than  a  certain 
value.^  The  volume  occupied  in  image  space  by  all 
distinct  image  points  within  a  radius,  cl,  divided  by 
the  volume,  V=Lj:LyL:,  occupied  by  each  image  point, 
gives  for  X p  in  this  case: 

Np'=  —(cty+—{dy+  —id),         (14) 

6V  32V  SV 

where : 

S=2{Lyy+LyL,-\-L,L,) 
LT  =  ^L,+  Ly+L,). 

Similar  arguments  for  the  cases  when  both  source  and 
receiver  are  in  a  comer,  at  the  center  of  a  wall,  or  at 
the  center  of  an  edge  show  that  all  these  cases  (except 
Case  1,  of  course)  can  be  expressed  by  the  formula : 


N 


1       [Tr{cty     Tr 

-+-iVxL.+  r,yLy+r,.U{dy 


-\-4{ViVvLxLy+r)y7j!LyL:-\-r]ir]j:LiLx)(d)    ,      (15) 


I 


where  tj^,  rjy,  and  ?)..  are  to  be  chosen  as  follows: 

r]x=  r]y=  7]t=  2,  source  and  receiver  in  comer ; 

»?i=T7y=2,  source  and   receiver  at   center  of  a 

r;t=  1  z-edge,  etc.; 

r]z=2        ,  source   and    receiver   at   center   of   a 

T;„=77i=l  y2-wall,  etc.; 

i]z=Vi/=Vz=  li  source  and  receiver  at  center  of  room. 

It  is  clear  that  for  a  given  room  the  formula  for  .Vp'  will 
depend  upon  the  position  of  both  the  source  and  re- 
ceiver. Further,  for  pulses  of  a  finite  length,  for  actual 
sources  and  receivers,  and  for  most  positions  of  source 
and  receiver,  exact  coincidences  will  be  rare,  and 
"almost  coincidences,"  with  interference  effects  between 
the  various  pulses,  will  be  the  rule  rather  than  the 
exception. 

3.  Pulse  Spacing  Statistics 

The  close  similarity  between  the  normal  frequency 
lattice  and  the  three-dimensional  image  source  lattice 
has  already  been  noted.  The  average  statistical  proper- 
ties of  pulses  bear  a  close  resemblance  to  the  equivalent 
properties  of  the  normal  frequencies.  Therefore,  simi- 
larities in  the  fluctuation  statistical  properties  common 
to  both  normal  frequencies  and  pulses  might  be 
expected. 

The  frequency  spacing  index  \p  has  been  defined  as 
the  mean  squared  ratio  of  actual  to  average  normal 
frequency  spaces,  for  a  specified  frequency  interval. 
This  index  has  been  evaluated  for  a  certain  class  of 
rectangular  rooms,  and  it  can  be  calculated  if  normal 
frequency  values  are  available.^  We  shall  see  next  that 
the  exact  analog  of  xp  for  pulse  spacing  (defined  as  the 
mean  squared  ratio  of  actual  to  average  time  intervals 
for  returning  pulses)  can  be  calculated  and  evaluated 
in  every  case  for  which  the  frequency  spacing  index 
index  can  be  calculated  or  evaluated,  provided  that 
source  and  receiver  are  maintained  in  one  corner. 

We  start  with  the  expression  for  the  dimensionless 
normal  frequency  of  the  /,  m,  nth  mode:^ 


I,  m,  11  =  0,  1,2-  ■  ■. 


(16) 


We  obtain  the  time  of  arrival  of  distinct  pulses  from 
the  image  sources  by  using  Eq.  (15)  for  the  case  where 
source  and  receiver  are  both  in  a  corner.  The  resulting 
equation  for  time  of  arrival  can  be  put  into  the  dimen- 
sionless form, 

r/4=hipq)-^[.P+impy+(nqyy,  (17) 

/,  w,  w  =  0,  1,  2-  •  •, 

where  t  is  as  defined  in  Eq.  (12). 

Notice  that  Eq.  (17),  which  applies  only  to  the  case 
where  source  and  receiver  lie  in  the  same  comer,  differs 
from  Eq.  (16)  for  nimn  only  through  the  appearance  of 
{p)~^  and  {q)~~^  in  place  of  p  and  q.  Since  p  and  q  enter 

'  R.  H.  Bolt,  "Normal  frequency  spacing  statistics,"  J.  Acous. 
Soc.  Am.  19,  79  (1947). 
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symmetrically  in  both  Eqs.  (16)  and  (17),  they  can  be 
interchanged.  Thus  m  order  to  convert  any  equation 
giving  the  frequency  spacing  index  into  a  form  which 
will  yield  the  pulse  spacing  index,  one  need  only  substi- 
tute p~^  for  q,  q"^  for  p,  and  r/4  for  ju.  Of  course,  if  the 
room  proportions  are  such  that  the  required  weighting 
factors  cannot  be  evaluated,  but  actual  pulse  spacing 
values  are  available,  for  example  from  experiment,  ip 
for  pulse  spacing  is  easily  evaluated  from  its  defining 
equation : 

1  B     /8Timn\- 

^AB  = E   (  — —  I    (5r>Av,  (18) 

V(6t)av/ 


TB—TA 


where  (5t)av/4  is  obtained  by  making  the  substitutions 
p-^{q)~^,  q-^ip)~^,  M~*i"/4  in  the  equation  giving  the 
average  spacing  between  adjacent  normal  frequencies.' 

4.  A  Derivation  of  the  Mean  Free  Path 

Pulse  analysis  yields  a  straightforward  derivation  of 
the  classical  mean  free  path.  We  assume  that  the  order 
numbers  are  so  large  that  the  coordinates  of  the  source 
and  receiver  can  be  neglected  (that  is,  source  and  re- 
ceiver can  be  considered  to  be  at  the  origin).  Using 
polar  coordinates  (r,  6,  </>)  we  can  express  the  number  of 
reflections  associated  with  each  cell  as : 


Nu 


cos<j>  cos6    sin0     cos</)  sinS 


[cos<^  COS0     sm 
1 


(19) 


We  next  obtain  an  average  value  of  the  number  of 
reflections  out  to  a  given  radius  r  =  ct,  by  averaging 
Eq.  (19)  over  all  angles: 
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Nlmnda 


N,     =- 

■*■  '  Imn 


So 


rS 


(20) 


where  da  =  r-  cos(f)dct)dd ,  So=  {'iirr'')/8,  S=2{LyLz-\-LzLi 
-\-LxLy),  and  the  integration  is  taken  over  one  octant 
of  space.  By  definition,  the  mean  free  path  is  equal  to 
the  total  distance  traveled  by  an  average  pulse  in  a 
given  time  /,  divided  by  the  number  of  reflections  of 
that  pulse  during  the  same  time.  Therefore : 


m.i.^.  =  r/Nimn  =  ^V/S. 
5.  Energy  in  the  Pulse 


(21) 


We  now  consider  the  energy  contained  in  a  sound 
pulse  and  follow  the  course  of  that  energy  as  the  sound 
becomes  dispersed  throughout  the  room  and  absorbed 
at  its  boundaries.  The  total  power  radiated  by  a  simple 
source  is  :^ 

wpv-Qo^    4t 

n  = =  -({p^W)  ergs/sec,  (22) 

2c         pc 

'  P.  M.  Morse,  Vibration  and  Sound  (McGraw-Hill  Book  Com- 
pany, Inc.,  New  York,  1948).  second  edition. 


where  Q^,  is  the  source  strength,  and  (/>o^)  is  the  mean 
square  sound  pressure  at  a  distance  r^  from  the  source 
in  a  free  field.  This  equation  applies  to  steady-state 
radiation.  It  is  also  valid  for  a  single  pulse  wave  train 
which  has  a  sufficiently  narrow  spectral  distribution. 
This  requirement  is  fairly  well  satisfied,  for  example, 
if  there  are  ten  or  more  waves  in  the  train,  and  if  the 
wave  amplitude  is  fairly  constant  throughout  the 
duration  of  the  pulse.  We  designate  the  length  of  this 
pulse  Tp  and  write  the  total  energy  contained  in  the 
pulse  as : 


■■  nTp=  (4Tr/pc)i{payo^)Tp  ergs. 


(23) 


The  energy  density  in  the  pulse  is  continually  diminish- 
ing as  the  pulse  radiates  outward.  At  any  instant  of 
time,  /,  after  emission,  the  volume  of  space  containing 
the  pulse  is : 


Vp=  (47r/3)c^[^-  {t-  rp)3]~47rc3/2rp. 


(24) 


We  next  consider  the  multiplicity  of  pulses  arriving  at 
the  receiver  from  all  directions  as  time  progresses  (still 
assuming  that  the  room  is  lossless).  We  take  the  number 
of  pulses  arriving  per  second,  as  given  by  Eq.  (11),  and 
multiply  this  by  the  energy  density  during  the  duration 
of  a  pulse  passage,  obtaining  the  energy  density  per 
second  (provided  the  pulses  add  incoherently) : 

£p  dNp     47r  iPo'W 

Wp= = ergs/cc-sec.        (25) 

Vp    di      pc      V 

If  we  now  multiply  \Vp  by  the  fraction  of  a  second 
occupied  by  the  individual  pulse,  we  obtain  the  average 
energy  density  in  a  room  of  volume  V : 


47r  {po'W 

Wr  = Tp  ergs/cc. 

pc       V 


(26) 


A  more  convenient  quantity  experimentally  is  the  mean 
square  pressure  in  the  room : 


(Pi^)  =  PcWr  = 


^Trc{{po^)ra^)Tp 
V 


(27) 


Suppose  we  take  a  pulse  of  rp=  20  msec,  duration 
which  generates  a  mean  square  pressure  of  one  dyne 
(74  db  pressure  level)  at  a  distance  of  one  meter  from 
the  source  in  a  free  field.  If  this  pulse  is  emitted  in  a 
room  of  10,000  cu.  ft.  volume,  we  find  that  the  average 
mean  square  pressure  throughout  the  room  after  the 
pulse  is  dispersed,  as  calculated  from  Eq.  (27),  is  about 
69  db,  or  only  5  db  less  than  the  sound  pressure  level 
in  the  original  pulse  as  measured  at  one  meter.  In 
practice,  the  sound  is  rapidly  dissipated  by  absorption 
and  may  be  canceled  out  by  interference  effects. 
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6.  A  Derivation  of  the  Eyring 
Reverberation  Equation 

Let  us  assume  thai  all  of  the  boundaries  of  the  room 
are  equally  absorptive,  and  that  their  absorption  is  less 
than  20  percent,  so  that  the  specular  reflection  ap- 
proximation is  valid.  Further,  let  us  suppose  that  the 
energy  in  a  pulse  is  diminished  upon  reflection  by  a 
factor  (1  — a),  where  a  is  ihe  average  over  all  angles  of 
the  free  wave  absorption  coefficient.  Then  the  pulse 
energy  associated  with  the  A'th  image  is: 

Ej,s  =  E,{l-a)^=Ep{l-ay^'S''^\  (28) 

which  leads  at  once  to  the  Eyring  reverberation  equa- 
tion: 

r6o  =  A'F/51n[l,''(l-a)]. 

This  is  the  average  long  term  reverberation  time  of  the 
room.  To  study  the  short  term  response  for  a  fairly 
hard-walled  room,  one  could  utilize  the  specular  reflec- 
tion approximation^' •*  for  each  of  the  first  few  reflections 
and  thus  obtain  in  detail  the  approximate  short  term 
response  of  the  room  to  a  particular  transient.  Further 
refinements  of  the  long  term  response  might  be  obtained 
by  treating  groups  of  images  statistically,  just  as 
modified  decay  equations  have  been  obtained  by  group- 
ing of  normal  modes.^ 

PULSE  MEASUREMENTS  IN  A  HARD-WALLED 
RECTANGULAR  ROOM 

1.  Experimental  Procedure 

In  order  to  obtain  some  simple  results  which  could  be 
interpreted  from  a  pulse  statistics  point  of  view,  short 
sound  pulses  were  produced  experimentally  in  a  hard- 
plaster-walled  rectangular  room,  and  photographs  of  the 


sound  arriving  at  a  microphone  were  made  using  a 
cathode-ray  oscillograph.  The  source  used  was  a  W.E. 
713A  receiver  unit  feeding  into  a  f-in.  diameter  brass 
tube  8  in.  long  packed  with  steel  wool  to  present  a  high 
acoustic  impedance  to  the  diaphragm.  The  effective 
source  was  the  end  of  the  tube,  which  was  small  and 
could  be  considered  approximately  a  simple  source. 
The  microphone  was  a  W.E.  633A  dynamic  type,  the 
output  of  which  was  passed  through  an  ERPI  RA-363 
octave  filter  and  through  the  amplifier  section  of  a 
ERIT  RA  277-F  sound  analyzer. 

The  output  of  the  analyzer  was  attenuated  logarith- 
mically by  a  Kay  Labs  Type  510-A  Logaten  and  then 
put  across  the  vertical  deflection  plates  of  a  DuMont 
247  cathode-ray  oscillograph.  A  pulsed  carrier  of  about 
2  msec,  duration  was  produced  by  mechanical  switching 
of  the  speaker  input  which  consisted  of  a  3600-c.p.s. 
signal  from  a  Hewlett  Packard  200  D  audio  oscillator 
amplified  through  a  Fairchild  audio  amplifier.  A  slow 
(approximately  42  cm/sec.  on  the  screen)  external  single 
sweep  on  the  oscillograph  was  activated  mechanically 
shortly  before  the  beginning  of  the  pulse.  Thus  the 
logarithm  of  the  acoustic  signal,  picked  up  by  the  micro- 
phone during  the  0.25-sec.  interval  after  the  pulse  was 
emitted,  was  recorded  linearly  on  the  oscillograph 
screen  as  a  linear  function  of  time,  and  was  photo- 
graphed. 

2.  Calibration  and  Auxiliary  Data 

The  horizontal  sweep  speed  was  calibrated  by  a 
30-c.p.s.  signal  direct  to  the  vertical  plates  from  the 
oscillator.  The  vertical  deflection  was  calibrated  in  5-db 
steps,  using  the  attenuator  pad  on  the  ERPI  analyzer. 
The  horizontal  sweep  was  found  to  be  linear  within 
experimental  error,  while  the  vertical  deflection  was 
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found  to  be  linear  in  db  with  a  maximum  deviation  of 
±1  db  over  the  main  working  range  of  50  db.  With 
somewhat  larger  variations  from  linearity  at  low  ampli- 
tudes, the  vertical  deflection  is  approximately  linear 
over  a  range  of  70  db,  which  includes  background  level 
in  the  pictures. 

The  following  assumptions  and  conventions  were 
adopted  in  interpretation  of  the  data:  (a)  All  pictures 
were  cahbrated  in  db  vs.  ft.,  since  time  in  the  room 
corresponds  to  distance  in  image  space;  (b)  it  was  as- 
sumed that  the  acoustic  pulse  emitted  by  the  source  in 
all  pictures  was  the  same  in  both  shape  and  height;  (c) 
a  decibel  reference  level  was  established  in  terms  of  the 
peak  amplitude  and  was  taken  as  zero  db  at  0.1  ft. 
from  the  source,  all  other  levels  thus  coming  out  nega- 
tive; (d)  the  velocity  of  sound  was  taken  as  1120  ft./sec. 

In  order  to  get  a  rough  check  on  the  validity  of  the 
assumption  that  the  actual  source  was  a  simple  source, 
two  sets  of  pictures  were  taken  in  which  source  and 
receiver  were  placed  near  the  center  of  the  room  (so  that 
direct  pulses  would  arrive  without  interference  from 
images)  and  were  spaced  several  different  distances 
apart  from  0.50  to  4.0  ft.  The  peak  levels  of  these  direct 
pulses  plotted  against  distance  between  source  and 
receiver  on  semi-log  paper  are  shown  in  Fig.  4.  The  best 
straight  line  fit  to  the  experimental  points  is  a  line  with 
inverse  square  slope.  This  in  turn  indicates  that  in- 
terpreting the  actual  source  as  a  simple  point  source  is  a 
fair  approximation.  Another  check  on  this  approxima- 
tion was  made  by  measuring  peak  pulse  amplitude  in 
db  at  a  constant  radius  from  the  source,  but  at  six 
different  angles  with  respect  to  the  tube,  from  0°  to 
180°.  These  measurements  showed  a  maximum  varia- 
tion of  ±2  db  in  the  peak  values. 

3.  Pulse  Shape  and  Energy 

Pulse  length  and  carrier  frequency  were  selected 
empirically.  The  pulse  length  was  chosen  short  enough 
compared  to  room  dimensions  to  give  at  least  several 
clearly  separated  "echoes,"  but  long  enough  to  include 
at  least  five  or  ten  cycles  of  the  carrier.  In  turn  the 
carrier   frequency    was   selected   as   high    as   possible 
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Fig.  5.  Photograph  for  a  non-degenerate  case  in  which  the 
distance  between  source  and  receiver  is  1'.  This  photograph  shows 
a  typical  direct  pulse  without  interference  from  image  pulses. 


Fig.  6.  Graphical  enlargement  of  the  direct  pulse  shown  in  Fig.  5, 
illustrating  the  general  shape  of  the  envelope. 

(without  encountering  large  air  absorption)  in  order 
that  the  speaker  would  operate  in  a  region  of  high 
efficiency  and  good  transient  response.  The  sweep  speed 
and  carrier  frequency  were  such  that  carrier  details 
are  barely  unresolvable,  although  individual  oscillations 
of  the  carrier  can  be  distinguished  in  some  places  in  the 
pictures. 

The  exact  nature  and  shape  of  the  pulse  was  subse- 
quently determined  by  analysis  of  the  pictures  (e.g.. 
Fig.  5).  The  average  pulse  dimensions  as  determined 
from  the  photographs  are  illustrated  in  Fig.  6,  which  is  a 
graphical  enlargement  of  tlie  direct  pulse  shown  in 
Fig.  5.  The  detailed  shape  of  the  pulse  is  represented  by 
a  stepped  curve  in  which  the  width  of  most  of  the  steps 
is  equivalent  to  several  cycles  of  the  carrier.  At  the 
onset  of  the  pulse  the  first  two  or  three  swings  of  the 
alternating  and  e.xponentially  increasing  carrier  are  just 
discernible  and  the  pulse  builds  up  to  maximum  ampli- 
tude in  about  two  cycles  of  the  3600-c.p.s.  carrier.  The 
pulse  is  almost  symmetric  but  its  peak  ampHtude  in  the 
positive  direction  is  measurably  greater  than  in  the 
negative  direction  for  this  particular  picture.  (Sym- 
metry of  the  pulse,  of  course,  depends  upon  the  phase 
of  the  carrier  at  the  instants  of  switching  it  on  and  off.) 
Following  the  peak,  the  pulse  decays  in  a  rigorousl)' 
exponential  manner  within  experimental  error  as  indi- 
cated in  Fig.  6  by  the  two  straight  lines  drawn  through 
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the  stepped  curves.  'J'he  pulse  decay  rate  is  about 
2.5X10^  db/sec.,  and  depends  principally  upon  the 
loudspeaker  characteristics.  The  background  noise 
showinji;  on  Fig.  5  is  mainly  a.c.  power  ripple.  Maximum 
background  level  for  all  pictures  was  —65  db. 

The  energy  in  a  pulsed  carrier  having  an  envelope 
that  rises  instantaneously  and  then  decays  exponen- 
tially is 

TTpu-ru-    W-+  27- 

Ep  =  — ergs,  (29) 

pc      7(w^+7^) 

where  />o  is  the  peak  pressure  in  dynes  at  ro  cm,  co  is  the 
angular  frequency,  and  7  is  the  exj)onential  decay  con- 


stant. For  the  pulse  sliown  here  7  =  580  and  co=72007r 
so  that  {y/u)y<^l.  Hence  Eq.  (29)  reduces  to 


E..= 


irPoW 


(30) 


pcy 


If  we  let  p,  be  the  r.m.s.  pressure  in  the  room  a  rela- 
tively long  time  after  emission  of  the  pulse  (i.e.,  when 
the  sound  is  fairly  uniformly  distributed),  then  the 
ratio  Ps/pa  is,  using  Eq.  (27) : 
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Figs.  7-12.  Series  of  pulse  photographs  showing  response  in  a  hard-walled  rectangular  room  when  the  source  is  placed  in 
a  corner  and  the  receiver  is  located  at  various  points  around  the  room.  Note  the  extreme  number  of  coincidences  when  the 
receiver  is  also  in  a  corner  (e.g.,  Figs.  8  and  11). 
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Fig.  13.  The  total 
number  of  distinct  pulses 
Np  arriving  at  the  re- 
ceiver, plotted  as  a 
function  of  time,  when 
the  source  and  receiver 
are  both  in  the  corner 
of  the  room,  as  in  Fig.  1 1 . 
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or,  using  Eq.  (30)  and  letting  Tp— §7: 


ps      /2Trro^CTp\' 
p,     \       V      ) 


(31) 


This  expression  neglects  dissipation  of  energy  in  the 
room  which  is  easily  taken  into  account  if  the  reverbera- 
tion time  of  the  room  is  known. 

4.  The  Period  of  Resolution 

A  period  of  resolution  tr  can  be  defined  as  the  time  at 
which  the  expected  interval  between  successive  pulses  is 
just  equal  to  the  effective  pulse  duration  Tp.  At  a  time 
tr  after  the  first  pulse,  one  would  no  longer  expect  to 
"see"  the  individual  echoes  in  the  clear,  but  rather  a 
smear  characterized  by  an  envelope  above  which  the 
peaks  of  the  pulses  would  occasionally  appear. 

An  expression  for  the  period  of  resolution  tr  in  terms 
of  the  pulse  width  r-p  can  be  obtained  from  the  correct 
(i.e.,  appropriate  for  the  existing  number  of  degenera- 
cies) expression  for  Np  or  N p  given  in  Eqs.  (10),  (13)- 
(15)  by  solving  the  difference  equation  for  /, : 


Source  and  receiver  in 
arbitrary  positions: 

Source  in  corner,  receiver 
in  arbitrary  position ; 

Source  and  receiver  in 
corners : 


/,  =  9.98X  10-3  sec. 


/r  =  26.3X10-3  sec. 


/r=  56.8X10-3  sec. 


Np{tr)-Np{tr-Tp)=\. 


(32) 


This  has  been  done  for  the  conditions  of  the  experi- 
ments performed  and  leads  to  the  following  values  for 
the  period  of  resolution: 


These  values  correspond  to  an  effective  pulse  duration  of 
Tj,=  1.74X10~^  sec,  which  is  the  pulse  width  at  an 
amplitude  equal  to  \/e  times  the  peak  amplitude.  The 
resolving  time  as  defined  is  depicted  in  each  of  the 
pulse  pictures  by  a  vertical  line  appearing,  in  every  case 
but  one,  to  the  right  of  the  initial  pulse.  The  individual 
pictures  will  be  discussed  in  detail  later. 

5.  Envelope  and  Decay  of  Unresolved  Sound 

The  height  of  the  vertical  fine  that  indicates  the 
resolving  time  has  been  adjusted  to  equal  the  r.m.s. 
sound  level  relative  to  the  peak  pulse  level  as  obtained 
from  Eq.  (31).  The  subsequent  decay  of  this  more  or 
less  continuous  sound  is  portrayed  by  the  two  slightly 
sloping  horizontal  lines  that  form  the  envelope  of  the 
unresolved  sound.  The  slopes  of  the  envelope  were  de- 
termined from  the  measured  reverberation  time  of  the 
room  at  3600  c.p.s. 

It  will  be  noted  that  the  size  of  the  envelope  is,  in  a 
number  of  cases,  considerably  too  large  (e.g..  Figs.  8, 
10,  11).  In  every  case  where  the  theoretical  envelope 
does  not  fit  the  photograph,  the  source,  or  both  the 
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THEORETICAL   AVERAGE    PULSE  RECEPTION 
RATE   PER   SECOND 

NUMBER   OF  PULSES    ACTUALLY    RECEIVED  IN 
INTERVAL   OF  154   MILLISECONDS 


Fig.  14.  The  aver- 
age and  actual  rates 
of  arrival  of  pulses  at 
the  receiver  when 
l)oth  source  and  re- 
ceiver are  in  the 
corner  of  the  room 
as  in  Fig.  11. 
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Figs.  15  and  16.  Pulse  responses  for  the  receiver  in  other  locations. 


source  and  receiver,  was  in  a  comer  position.  Since  the 
source  was  of  finite  dimensions  it  was  never  exactly  in  a 
comer  but  usually  from  0.2  to  0.5  of  a  wave-length  from 
one  or  more  of  the  walls.  Hence  there  was  interference 
between  the  source  and  its  images,  resulting  in  a  strongly 
directional  elective  source,  and,  consequently,  in  direc- 
tional image  clusters.  Thus  for  the  source  in  a  corner 
and  a  particular  location  of  the  receiver,  certain  image 
sources  will  contribute  ver>-  little  efTective  energy  to 
the  receiver.  This  amounts  to  saying  that  the  source 
effectively  puts  less  energy  into  the  room.  This  point 
will  be  made  more  evident  from  the  discussion  of  the 
photographs. 

6.  Discussion  of  Pulse  Photographs 

The  photographs  were  taken  in  a  room  of  dimensions 
Lx=23.0,  Ly=UA  and  Z,  =  8.44  ft.  To  specify  the 
location  of  the  source  and  receiver,  we  shall  use  (A',  F,  Z) 
for  the  source  position  and  {U,  V,  IF)  for  the  receiver 
position,  distances  being  measured  in  feet. 

In  Fig.  5,  the  source  coordinates  are  (11,  7.2,  3.9  ft.) 
and  the  receiver  was  located  at  (12,  7.2,  3.9  ft.).  This 
position  is  non-degenerate  and  so  the  resolving  time 


shown  is  lr=9.98  msec.  The  envelope  is  seen  to  be 
slightly  too  large.  This  can  perhaps  be  accounted  for 
by  the  fact  that  both  receiver  and  source  are  near  the 
center  of  the  room  so  that  images  from  opposite  pairs 
of  walls  interfere  with  each  other  consistently. 

For  Fig.  7,  (XYZ)  are  (1,  1.5,  0.5  in.)  and  {UVW)  are 
(5.75,  3.35,  2.11  ft.).  The  source  is  thus  effectively 
eightfold  degenerate,  being  in  the  comer,  so  the  resolv- 
ing time  is  tr=  26.3  msec.  It  will  be  noted  that  the  direct 
pulse  is  exceeded  in  height  by  two  succeeding  pulses. 
The  peak  amplitudes  of  the  pulses  indicated  by  the 
squares  are  plotted  on  Fig.  4.  This  sequence  corresponds 
to  pulses  from  the  successive  images  at  (0, 0, 0), 
(0,  0,  2L,),  (0,  0,  4Z-,),  and  (0,  0,  6L,),  each  of  which  is 
a  cluster  of  eight  image  sources.  The  differences  in 
amphtude  show  clearly  that  each  image  cluster  has 
distinct  directionality,  with  especially  great  variations 
in  the  A',  V  plane.  The  long-dashed  line  on  Fig.  4  is  the 
peak  level  one  would  expect  for  incoherent  addition  of 
16  sources  (i.e.,  12  db  above  peak  for  a  single  source). 
It  is  seen  that  the  direct  sound  is  12  db  below  this  line. 
However,  the  sound  from  the  image  sources  hes  along 
a  line  14  db  above  the  peak  for  a  single  source,  which  is, 


339 


PULSE     STATISTICS 


in  turn,  4  db  below  the  level  expected  if  the  eight  sources 
in  each  cluster  all  added  coherently. 

Figures  8-12  are  a  series  in  which  (XYZ)  are  still 
(1,  1.5,  0.5  in.)  and  in  which  U=l  in.,  W=l  in.  for 
all  photographs.  In  Fig.  11,  V  =  Ly—\;  in  Fig.  12, 
V  =  3Ly/4:;  in  Fig.  13,  V  —  Ly/2;  and  in  Figs.  14  and  15, 
V=l  in.  The  directionahty  of  the  image  clusters  is 
clearly  evident  in  Fig.  8.  Here  there  is  no  distinct  "di- 
rect" pulse  since  the  pulse  from  the  image  at  (0,  2Ly,  0) 
arrives  simultaneously  with  the  pulse  from  the  actual 
source.  Evidently  these  two  pulses  interfere  since  the 
height  of  the  first  pulse  is  small.  The  series  of  tall 
pulses  is  again  from  the  images  at  (0,  0,  2Z,;),  (0,  0,  4ZJ 
etc.,  each  of  which  arrives  simultaneously  with  corre- 
sponding pulses  from  the  images  at  (0,  2Ly,  IL^), 
(0,  2Ly,  4L2),  and  their  reflections  in  the  X,  Y  plane. 
Each  of  these  pulses  then  represents  the  simultaneous 
contributions  of  32  image  sources.  As  is  evident  by 
comparison  of  the  peak  heights  with  those  in  Fig.  7, 
the  sources  do  not  add  coherently,  however.  The 
marked  effect  of  interference  is  thus  responsible  for  the 
low  value  of  the  actual  envelope  in  this  picture.  The 
symmetry  of  the  images  corresponds  closely  to  the  sym- 
metry when  both  source  and  microphone  are  in  the 
same  corner  so  the  resolving  time  indicated  is  /r=56.8 
msec. 

Figures  9  and  10  show  effects  similar  to  those  just 
discussed  for  the  other  positions  of  the  microphone 
mentioned  previously.  The  resolving  time  for  these  is 
again  /r=26.3  msec,  since  the  microphone  position  is  no 
longer  "degenerate."  The  influence  of  image-source 
directionality  continues  to  be  apparent. 

In  Figs.  11  and  12  both  source  and  receiver  are  in  the 
corner.  Figure  12  serves  to  indicate  the  reproducibiUty 
of  Fig.  11  which  was  possible  with  the  experimental 
arrangements  employed.  The  resolving  time  here  is 
ir=  56.8  msec.  The  triangles  refer  to  the  points  shown  on 
Fig.  4.  These  pulses  each  represent  the  simultaneous 
contributions  of  the  images  at  (0,  0,  ±2Lj),  (0,  0,  ±4^2) 
etc.  As  is  evident  from  Fig.  4,  the  peaks  very  nearly 
follow  the  inverse  square  law.  The  peak  amplitude, 
however,  corresponds  to  9  db  above  peak  for  a  single 
source,  whereas  each  pulse  comes  from  16  sources. 
The  resulting  peak  in  this  case  is  less  than  would  be 
expected  for  incoherent  addition.  Consequently  it  is 
again  not  surprising  that  the  actual  envelope  is  below 
the  predicted  value.  The  lines  drawn  on  the  lower  half 
of  the  photograph  indicate  the  calculated  times  of  ar- 
rival of  distinct  pulses.  Figure  13  shows  the  total  num- 
ber of  distinct  pulses  Np  arriving  up  to  a  time  /,  as 
calculated  from  Eq.  (15)  for  this  location  of  source  and 
receiver;  and  Fig.  14  shows  the  average  rate  of  arrival, 
and  a  graphical  portrayal  of  the  actual  pulses  arriving. 
Each  line  on  Fig.  11  corresponds  to  a  pulse  on  Fig.  14. 
It  is  seen  that  agreement  is  quite  good. 

In  Fig.  15  the  source  is  still  in  the  comer  but  the 
receiver  is  at  (|Lx,  jLy,  1  in.).  This  picture  should  be 
compared  with  Fig.  10. 


In  Fig.  16  the  source  remains  at  (1,  1.5,  0.5  in.)  while 
the  receiver  is  placed  in  the  opposite  corner  (nearly)  at 
(Lx,  Ly,  0).  Here  the  pulse  from  the  source  arrives  al- 
most simultaneously  with  the  pulses  from  the  three 
image  clusters  at  (0,  2L„,  0)  (2Lj.,  2Ly,  0),  and  (2Lx,  0, 0). 
The  resulting  interference  is  clearly  evident  from  the 
appearance  of  the  first  pulse.  In  this  case  the  symmetry 
is  again  almost  the  same  as  that  for  both  source  and 
receiver  in  the  same  comer.  Therefore  /r=56.8  msec. 
is  used.  Agreement  of  actual  and  calculated  envelopes 
is  seen  to  be  fairly  good. 

Although  the  present  experiments  are  restricted  to  an 
idealized  room,  a  close  connection  can  already  be  seen 
between  these  results  in  a  small  room  and  pulse  studies 
in  large  auditoriums.  For  example,  in  Fig.  17  the  source 
is  at  (1.5,  6.7,  1  ft.)  and  the  receiver  at  (fZ-x,  jL„,  9  in.). 
These  locations  correspond  to  typical  locations  of  a 
speaker  and  listener  in  a  large  rectangular  hall.  If  we 
consider  a  scaling  factor  of  about  5,  the  dimensions  of 
this  hall  would  be  115X67X42.2  ft.,  and  our  pulse 
would  have  a  corresponding  width  (at  1/e  peak  ampH- 
tude)  of  8  msec.  The  scaled  carrier  frequency  would  be 
about  700  c.p.s.  The  resolving  time  is  seen  to  be  shorter 
than  the  time  actually  required  for  the  direct  pulse  to 
reach  the  receiver,  and  the  character  of  the  response 
shows  that  a  "smear"  sets  in  immediately.  The  actual 
and  calculated  envelopes  are  seen  to  agree  well  in  this 
case,  since  no  degeneracies  are  present. 

Figures  18  and  19  can  be  compared  qualitatively  with 
Fig.  17.  These  were  obtained  in  a  motion  picture  theater 
of  about  250,000  cu.  ft.  volume.  Pulses  were  produced 
by  applying  a  short  pulsed  carrier  signal  to  the  theater 
loudspeaker.  In  Fig.  18  the  microphone  was  at  a  seat 
in  the  side  section  of  the  main  floor  about  two-thirds  of 
the  way  back  from  the  stage,  under  the  balcony.  The 
pulse  length  was  5  msec,  and  the  carrier  frequency 
2800  c.p.s. 

In  Fig.  19  the  microphone  was  in  the  center  section 
of  the  main  floor,  well  in  front  of  the  balcony  and  some- 
what off  the  center  line.  The  pulse  length  was  again  5 
msec,  and  the  carrier  frequency  was  1500  c.p.s.  Quahta- 
tive  observations  indicated  that  the  "hour  glass"  bulge 
in  Fig.  19  can  be  correlated  subjectively  with  a  "slap 
echo"  from  the  hard  but  irregular  rear  wall.  This  bulge 
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Fig.  17.  Pulse  response  for  non-degenerate  locations  of  source 
and  receiver  which  were  chosen  to  correspond  respectively  to  the 
positions  of  a  speaker  and  listener  in  a  large  rectangular  hall. 
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Figs.  18  and  19.  Tj-jDical  pulse  responses  in  a  motion  picture  theater  of  about  250,(XX)  cu.  ft.  volume. 
Note  qualitative  similarity  with  Fig.  17. 


was  characteristic  of  pulse  responses  at  this  location  for 
nearly  all  carrier  frequencies.  A  more  detailed  analysis 
of  these  and  similar  photographs  is  in  progress  in  an 
attempt  to  correlate  various  features  of  the  pulse 
response  with  the  results  of  subjective  listening  tests. 

CONCLUSIONS 

Most  of  the  theoretical  and  experimental  work  pre- 
sented in  this  paper  is  directly  appHcable  only  to  hard- 
walled  rectangular  rooms.  However,  as  has  been 
pointed  out,  appropriate  modifications  of  the  images 
may  be  made  when  the  walls  are  not  hard  and  in  some 
cases  modifications  for  the  first  few  images  do  not 
involve  too  great  computational  difficulties. 

It  appears  that  fairly  detailed  experimental  and  theo- 
retical investigations  of  the  short  term  transient  re- 


sponse (e.g.,  the  first  100-200  msec.)  for  rooms  having 
various  shapes  and  absorptive  treatments  should  prove 
highly  instructive.  At  the  same  time,  the  statistical 
methods  outlined  in  this  paper  can  be  applied  to  the 
analysis  of  the  fluctuation  characteristics  of  the  long 
term  response. 

The  exj)erimental  results  presented  show  that  ap- 
propriate transients  can  be  easily  produced  and  ob- 
served, and,  in  simple  cases  at  least,  correlated  with 
boundary  conditions.  Also  it  appears  that  pulse  methods 
are  readily  adaptable  to  scale  model  experiments. 
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